AD-A208  036 


@ 


ARL-FLIGHT-MECH-T.V,.403 


Af-{ -005-586 


DEPARTMENT  OF  DEFENCE 

DEFENCE  SCIENCE  AND  TECHNOLOGY  ORGANISATION 
AERONAUTICAL  RESEARCH  LABORATORY 

MELBOURNE,  VICTORIA 


Flight  Mechanics  Technical  Memorandum  403 

INCORPORATION  OF  VORTEX  LINE  AND  VORTEX  RING  HOVER 
WAKE  MODELS  INTO  A  COMPREHENSIVE  ROTORCRAFT 
ANALYSIS  CODE  (U) 


DTIC 

^  ‘-I  .ECTE  I 
,  MAY  2  4 1989 


R.  ToffoleMo,  N.E.  Gilbert, 
5.  Hill,  K.R.  Reddy 


Approved  for  Public  Release 


r;~.  ~r  7 -  ■ 

Atm. -■>'40  re  '.j  | 


(C)  COMMONWEALTH  OF  A'JSTRALIA  I985» 


8  O 


JANUARY  190? 


0  -  0 


This  work  is  copyright.  Apart  from  any  fair  dealing  for  the  purpose  of  study, 
research,  criticism  or  review,  as  permitted  under  the  Copyright  Act,  no  part 
may  be  reproduced  by  any  process  without  written  permission.  Copyright  is  the 
responsibility  of  the  Director  Publishing  and  Marketing,  AGPS.  Inquiries  should 
be  directed  to  the  Manager,  AGPS  Press,  Australian  Government  Publishing 
Service,  GPO  Box  84,  Canberra,  ACT  2601. 


AR-005-586 

DEPARTMENT  OF  DEFENCE 
DEFENCE  SCIENCE  AND  TECHNOLOGY  ORGANISATION 
AERONAUTICAL  RESEARCH  UBORATORY 


Flight  Mechanics  Technical  Memorandum  408 

INCORPORATION  OF  VORTEX  LINE  AND  VORTEX  RING  HOVER 
WAKE  MODELS  INTO  A  COMPREHENSIVE  ROTORCRAFT 
ANALYSIS  CODE 

by 

R.  TOFFOLETTO,  N.  E.  GILBERT,  S.  HILL,  and  K.  R.  REDDY 
SUMMARY 

The  incorporation  of  simplified  hover  wake  models  into  the  comprehensive  rotorcraft 
analysis  code  CAMRAD  is  described  and  examples  are  given  on  their  use.  The 
axisymmetric  models,  in  which  vortices  are  represented  by  either  straight  lines  or  rings,  are 
a  more  generalized  form  of  the  free  wake  models  of  R.  T.  Miller  at  MIT.  with  the  wake 
geometry  also  able  to  be  prescribed.  Incorporation  has  allowed  access  to  the  tabular 
representation  in  CAMRAD  of  airfoil  section  characteristics  as  functions  of  angle  of  attack 
and  Mach  number,  and  has  broadened  the  range  of  rotor  wake  models  in  the  code  to  include 
a  free  wake  hover  model  that  does  not  have  the  convergence  problems  of  the  existing  free 
wake  model  when  used  for  hover. 
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blade  chord  ** 

blade  lift  coefficient  ** 

ratio  of  thrust  coefficient  to  rotor  solidity 

core  burst  test  parameter 

factor  introducing  lag  in  solution 

empirical  scale  factor  for  concentrated  vortices  inboard  of  tip  vortex 

axial  settling  rates  of  tip  vonex  before  and  after  passage  of  following  blade 

radial  contraction  parameters  for  tip  vortex 

number  of  aerodynamic  segments 

blade  tip  Mach  number 

number  of  blades 

radial  and  axial  displacement  coordinates  (origin  at  rotor  hub  centre;  z  positive 
down) 

r  at  mid-points  of  aerodynamic  segments,  for  i  =  1 M 
r  at  edges  of  aerodynamic  segments  (from  root  to  iip),  for  i  =  1 M+1 
burst  vortex  core  radius 
vortex  core  radius 

vortex  core  radius  limited  to  a  minimum  of  0.005 

unburst  vortex  core  radius 

number  of  concentrated  vortices  along  the  blade 

number  of  vortex  line  or  ring  levels  in  intermediate  wake 

net  radial  and  axial  induced  velocity  components 

radial  and  axial  induced  velocity  components  due  to  far  wake 

radial  and  axial  induced  velocity  components  due  to  intermediate  wake 

downwash  at  blade  due  to  trailing  near  wake  ** 

self-induced  downwash  at  blade  ** 

blade  angle  of  attack  “ 

blade  bound  circulation  ** 

incremental  change  in  appropriate  quantity 

tolerance  for  induced  velocity  convergence 

blade  pitch  angle  ** 

longitudinal,  lateral,  and  vertical  induced  velocity  components  (funct's  of  r,  y) 
blade  inflow  angle  ** 
blade  azimuth  angle 


All  quantities  are  dimensionless  (based  on  density,  rotor  rotational  speed,  and  rotor  radius).  Quantities 
used  only  locally  to  simplify  expressions  are  not  included  here 
* '  Function  of  r 

f  Function  of  r,  z 


Subscripts 

m,n  as  for  subscipt  (s,t)  but  at  a  source  of  induced  velocity,  i.e.  when  calculating 

induced  velocity  at  (r,,, ,  z, ,),  contributions  are  summed  over  (m,n) 
max  maximum  value 

new  value  at  current  iteration 

old  value  at  previous  iteration 

s,i  value  at  concentrated  vortex  number  s  (from  tip),  where  1  <  $  <  S,  and  at  vortex 

line  or  ring  level  number  t  at  which  induced  velocity  is  to  be  calculated  (i.e. 
object).  Note:  1  <  s  <  S  and  1  <  t  <  T  in  intermediate  wake;  t  =  0  at  blade,  in 
which  case  subscript  is  dropped  (e.g.  r,o  s  r,);  t  =  T  +  1  for  far  wake 


1.  INTRODUCTION 


In  response  to  requests  from  the  Australian  Services  to  evaluate  performance  character¬ 
istics,  especially  for  hover,  of  helicopters  presently  operated,  as  well  as  those  being 
considered  for  procurement.  Aeronautical  Research  Laboratory  (ARL)  has  developed  an 
analysis  capability  in  the  area  of  hovering  rotor  aerodynamics  which  includes  both  inhouse 
and  acquired  codes . 

In  1987,  Reddy  and  Gilbert  compared  predicted  hover  performance  with  flight  data  for 
four  helicopters.*  Comparisons  were  also  made  of  main  rotor  blade  loading  distribution  for 
one  of  the  helicopters,  a  Sikorsky  S-58  (equivalent  to  Westland  Wessex).  Predictions  were 
obtained  using  three  nonuniform  inflow  rotor  wake  models  and  a  uniform  inflow  model 
based  on  momentum  theory.  The  nonuniform  wake  models  used  were  a 

•  helical  vortex  lattice  prescribed  wake  model, 

•  vortex  line  prescribed  wake  model,  and 

•  vortex  ring  free  wake  model. 

The  first  of  these  models  is  incorporated  in  CAMRAD  (Comprehensive  Analytical 
Model  of  Rotorcraft  Aerodynamics  and  Dynamics),  a  code  developed  by  Johnson  at  Ames 
Research  Center.^'^.'*'^  The  code,  which  was  acquired  in  1984  as  part  of  a  cooperative 
program  with  the  US  Army,  uses  straight-line  vortex  elements  joined  in  the  form  of  a 
helical  vortex  lattice  to  represent  the  trailed  and  shed  vorticity.  CAMRAD  was  also  used  to 
provide  uniform  inflow  predictions  (based  on  momentum  theory),  a  preliminary  process  in 
determining  the  trimmed  prescribed  wake  solution. 

The  second  model,  which  uses  infinite  and  semi-infinite  straight-line  vortices,  was 
developed  by  Reddy  *7  at  ARL  independently  of  very  similar  work  by  Miller  at 
Massachusetts  Institute  of  Technology  (MIT)  and  Beddoes  (unpublished)  at  Westland 
Helicopters  Limited  (WHL).  An  acquired  free  wake  hover  code  that  was  implemented  at 
ARL  by  Hill  and  Reddy(unpublished)  was  used  as  the  third  model.  The  method,  which 
generally  represents  infinite  line  vortices  by  rings,  is  based  on  one  of  the  variations  in 
Miller’s  method. 

To  investigate  the  consistency  of  these  methods,  the  same  parameters  and  empirical 
corrections  were  used  within  each  model  in  applying  to  each  helicopter,  and  where 
possible,  consistency  of  appropriate  quantities  between  models  was  also  maintained.  The 
major  identified  inconsistency  between  the  wake  models  was  the  different  representation  of 
airfoil  characteristics.  In  CAMRAD,  the  two-dimensional  airfoil  section  characteristics  are 
represented  in  tabular  form  as  functions  of  angle  of  attack  and  Mach  number,  compressibil¬ 
ity  effects  are  therefore  effectively  incorporated.  In  the  more  simplified  vortex  line  and 
vortex  ring  models,  the  characteristics  are  represented  by  a  constant  lift  curve  slope  and  a 
quadratic  drag  polar  without  corrections  for  compressibility.  It  was  planned  therefore  to 
incorporate  the  vortex  line  and  vortex  ring  models  into  CAMRAD,  principally  to  allow  the 
two-dimensional  airfoil  data  to  be  available  to  these  simpler  models. 

Since  Miller's  models  are  formulated  in  a  way  that  allows  a  free  wake  geometry  for 
both  the  vortex  line  and  vortex  ring  models,  it  was  decided  to  incorporate  his  models,  but  in 
a  more  generalized  form,  allowing  the  geometry  to  also  be  prescribed  using  the  options  in 
CAMRAD,  as  well  as  some  additional  features.  The  main  purpose  of  this  repon  is  to 
document  the  model  formulation  used  and  to  provide  the  information  necessary  to  run  the 
models. 
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2 .  WAKE  MODELS 

Comprehensive  codes  such  as  CAMRAD  allow  the  analysis  of  a  complete  rotorcraft, 
usually  allowing  for  two  separate  rotors.  However,  it  is  generally  assumed  sufficient  to 
consider  only  the  main  rotor  in  the  case  of  hover  for  a  conventional  helicopter  with  a  single 
main  rotor  and  anti-torque  tail  rotor.  For  performance  predictions,  standard  estimates  are 
then  made  for  the  power  requirements  of  the  tail  rotor,  accessories  and  transmission,  and 
aerodynamic  interference.  By  assuming  an  axisymmetric  wake,  the  harmonics  of  blade 
motion  and  the  shed  wake  can  be  neglected,  and  only  collective  control  needs  to  be  adjusted 
to  trim  to  a  specified  thrust. 

The  helical  vortex  lattice  model  in  CAMRAD  follows  the  common  approach  of  closely 
tracing  the  three-dimensional  helical  shape  of  both  the  strong  tip  vortex  and  inboard  vortex 
sheet.  Unfortunately,  this  apparently  straight-forward  approach  results  in  a  computation 
process  that  is  complex  and  computationally  demanding.  This  is  especially  so  for  the  hover 
case  where  there  is  no  large  uniform  relative  wind  due  to  the  translational  velocity  of  the 
helicopter.  This  means  that  the  wake  is  not  swept  away  from  the  rotor  and  hence  more  of 
the  wake  must  be  considered.  It  also  means  that  wake  induced  velocities  are  the  only 
velocities  present,  resulting  in  a  greater  sensitivity  of  the  wake  geometry  to  changes  in 
induced  velocity.  This  increased  sensitivity  leads  to  instabilities  and  slow  wake 
convergence  if  free  wake  models  are  used.'® 

The  simplified  axisymmetric  methods  described  here  are  an  attempt  to  overcome  the 
above  problems  for  the  hover  case.  The  basis  of  these  methods  is  that  the  continuously 
descending  helix  is  represented  by  vertically  separated  horizontal  vortex  lines  (for  a 
horizontal  rotor  disc),  which  are  either  straight  or  circular.  The  wake  is  divided  into  three 
regions,  which  are  defined  as  near,  intermediate,  and  far. 


Root  vortices  Tip  vortices 


a)  plan  view 


Rotor  blade 
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Far  wake 


b)  side  view 


Fig.  1  Vortex  Line  Wake  Model 
with  Concentrated  Far  Wake 


Fig.  2  Vortex  Rmg  Wake  Model 
with  Distributed  Far  Wake 
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The  three  wake  regions  are  illustrated  in  Figs  1  and  2  for  the  vortex  line  and  vortex  ring 
models  respectively.  These  are  first  described  in  relation  to  a  prescribed  wake  geometry 
using  either  of  the  options  in  CAMRAD,  i.e.  (a)  Landgrebe"  or  (b)  Kocurek  and 
Tangier.  *2  Both  are  based  on  model  rotor  flow  visualization  data. 

For  both  vortex  line  and  vortex  ring  models,  the  near  wake  is  represented  by  semi¬ 
infinite  straight  lines  attached  to,  and  in  the  plane  of,  the  blade,  with  a  greater  concentration 
towards  the  tip.  Based  on  the  observations  in  Kocurek  and  Tangier's  experiments  of  four 
well  defined  tip  vortices  below  the  blade,  the  intermediate  wake  is  represented  at  each  of  the 
corresponding  four  axial  levels  by  either  two  straight  infinite  vortex  lines  (Fig.  1)  or  two 
vortex  rings  (Fig.  2).  However,  the  generalized  manner  in  which  these  methods  are 
implemented  here  allows  the  number  of  these  levels  to  be  varied  (up  to  36).  One  of  these 
'concentrated'  rolled-up  vortices  is  located  at  the  outside  boundary  of  the  prescribed 
contracted  wake  and  represents  the  strong  tip  vortex,  and  the  other  is  located  directly 
beneath  the  root  cutout  and  represents  the  inboard  vorticity. 

Below  the  intermediate  wake  region,  it  is  believed  that  the  wake  is  unstable;  the  tip 
vortices  undergo  viscous  dissipation,  resulting  in  wake  expansion.  To  account  for  this 
region,  which  is  still  close  enough  to  the  rotor  disc  to  induce  significant  inflow,  Kocurek 
and  Tangier  proposed  a  vortex  ring  with  radius  equal  to  the  rotor  radius,  axial  location  at 
the  same  level  as  the  fourth  tip  vortex  beneath  the  rotor,  and  strength  of  four  times  that  of 
the  tip  vortex.  This  concept  is  adopted  as  an  option  for  the  far  wake  (referred  to  as 
'concentrated  far  wake')  in  the  form  of  either  a  ring  for  the  vortex  ring  model,  or  an  infinite 
straight  lir.e  replacing  the  ring  for  the  vortex  line  model,  as  shown  in  Fig.  1 .  The  other 
option  for  the  far  wake  included  is  the  one  given  by  Miller  (referred  to  as  'distributed', 
'sheet',  or  'distributed  sheet')  using  two  semi-infinite  vortex  planes  (for  vortex  line)  or 
cylinders  (for  vortex  ring  -  as  in  Fig.  2)  with  strength  determined  by  the  geometry  of  the 
intermediate  wake,  and  positioned  one  wake  spacing  below  each  of  the  last  inner  and  outer 
rings  of  the  intermediate  region.  This  latter  option  is  the  only  one  used  for  the  free  wake 
method.  For  the  prescribed  wake  method,  the  far  wake  may  be  neglected. 

When  the  wake  geometry  is  allowed  to  be  free  in  Miller's  simplified  models,  the 
difficulties  of  convergence  typical  of  vortex  lattice  models  are  not  generally  experienced. 
The  radius  and  axial  spacing  of  each  concentrated  vortex  in  the  intermediate  wake,  with  its 
consequent  effect  on  the  distributed  far  wake,  is  determined  by  the  requirements  for 
equilibrium  of  the  velocities,  this  being  the  essence  of  the  free  wake  method.  In  Miller's 
method,  up  to  three  concentrated  vortices  are  allowed  in  the  intermediate  wake  though  only 
two  are  used  in  the  prescribed  intermediate  wake  here  and  in  Ref.  1  by  Reddy  and  Gilbert 
when  using  Miller's  vortex  ring  free  wake  model.  The  generalized  formulation  here  allows 
this  number  to  be  increased  up  to  ten. 

The  computational  procedure  is  outlined  in  Appendix  A,  with  the  prescribed  wake 
method  incorporated  as  part  of  the  complete  method,  and  expressions  for  the  velocity 
components  induced  by  wake  vortex  elements  are  given  in  Appendix  B  (see  Ref.  8  for 
derivations).  Block  diagrams  showing  the  separate  structures  of  the  free  and  prescribed 
wake  methods  are  given  in  Appendices  C  and  D  respectively. 


3 


3.  PROGRAM  MODIFICATIONS 

Modifications  to  the  standard  VAX  780  version  of  CAMRAD  are  given.  Because  the 
new  axisymmetric  models  are  intended  to  be  used  for  a  single  rotor  configuration, 
modifications  made  to  Rotor  1  subprograms  are  not  similarly  made  to  Rotor  2 
subprograms. 

The  following  subprograms  (each  stored  as  a  separate  Fortran  file,  e.g. 
GEOMRl.FOR)  in  CAMRAD  (see  Ref.  5,  Part  II)  have  been  modified  (sec  Appendix  E): 


GEOMRl 

Calculate  wake  geometry  distortion 

INPTWl 

Read  wake  namelist 

PRNTWl 

Print  wake  input  data 

RAMF 

Calculate  rotor/airframe  periodic  motion  and  forces 

TRIM 

Trim 

TIMER 

F*rogram  timer 

Changes  to  the  VAX  VMS  operating  system  since  1984  have  resulted  in  the  output  of 
null  component  CPU  times  at  the  end  of  the  CAMRAD  output  file.  Modifications  to 
TIMER,  while  not  necessary  for  implementation  of  the  models,  are  therefore  included  in 
Appendix  E. 

The  following  added  subprograms  (all  included  in  the  file  WAKERl. FOR)  form  the 
basis  of  the  new  models  (see  Appendix  F): 

WAKERl 

Determine  induced  velocity  at  rotor  using  vortex  line  or  ring 
model 

VTXIF 

Calculate  induced  velocity  in  intermediate  wake  due  to 
intermediate  and  far  wake 

IVTERP 

Calculate  induced  velocity  along  blade  at  concentrated 
vortices 

ILINE 

Evaluate  expressions  for  velocity  induced  by  vortex  line  in 
intermediate  wake 

IRING 

Evaluate  elliptic  integral  expressions  for  velocity  induced  by 
vortex  ring  in  intermediate  wake 

ELLIPCON 

Calculate  constants  used  in  elliptic  integral  expressions 

FRING 

Evaluate  elliptic  integral  expressions  for  velocity  induced  by 
semi-infinite  vortex  cylinder  in  far  wake 

FLINE 

Evaluate  elliptic  integral  expressions  for  velocity  induced  by 
semi-infinite  vortex  sheet  in  far  wake 

PRESWG 

Determine  prescribed  wake  geometry 

The  above  modified  files  and  added  file  are  first  compiled.  After  then  obtaining  the 
main  program  object  file  CAMRAD.OBJ  and  library  object  file  CAMRAD.OLB  containing 
all  original  compiled  subprograms,  the  new  executable  file  CAMRAD.EXE  is  given  on 
typing 

SLINK  CAMRAD, GEOMRl, INPTWl, PRNTW1,BAMF, TRIM, TIMER, WAKERl, CAMRAD/LIB 


Nine  new  input  variables,  all  in  namelist  NLWAKE,  have  been  added  (see  Appendix  G 
for  description  and  default  values).  Also  included  in  Appendix  G  are  some  comments  on 
existing  CAMRAD  variables. 

4.  TEST  CASES 

To  demonstrate  the  various  model  options  and  illustrate  the  effect  of  including 
compressibility,  comparisons  are  made  of  blade  loading  distribution  for  the  S-58  using 
Scheiman’s  test  data'^  as  in  Ref.  1 .  Main  rotor  perfoimance  and  blade  loads  are  obtained 
from  CAMRAD  using  the  basic  H-34  (i.e.  S-58)  data  deck  and  NACA  0012  airfoil  tables. 

Fig.  3  shows  the  effect  of  compressibility  using  the  vortex  ring,  free  wake  model. 
Since  the  tip  Mach  number  is  relatively  low  (=  0.56),  the  effect  is  only  minimal  in  this 
example.  For  the  'compressible'  case  in  Fig.  3,  the  command  file  and  resulting  output  file 
(the  latter  in  abbreviated  form)  are  given  in  Appendices  H  and  J  respectively  .  The  vonex 
ring,  free  wake  model  is  selected  by  setting  OPMODL  =  2  and  LEVEL(l)  =  2.  Two 
rolled-up  vortices  (NIBV  =  2)  and  four  vortex  levels  in  the  intermediate  wake  (NIVL  =  4) 
are  specified.  The  free  wake  model  only  allows  the  distributed  sheet  far  wake  model.  The 
empirical  factor  scali.ig  the  rolled-up  concentrated  vortices  inboard  of  the  one  at  the  tip  is  set 
to  the  default  value  of  0.6.  By  setting  inputs  DEBUG(14)  and  DEBLIG(24)  to  1 ,  additional 
information  on  the  induced  velocity  and  free  wake  geometry  is  printed.  Because  the 
.solution  is  independent  of  azimuth,  the  number  of  azimuth  steps  per  revolution  (MPSl  and 
also  MPSIR)  is  set  to  the  minimum  value  of  4,  the  number  of  blades. 


Fig.  3  Effect  of  Compressibility  on  Blade  Load  Distribution  for  S-58  using 
Vortex  Ring,  Free  Wake  Model  (Ct/o  =  0.0817) 

Operating  conditions  and  main  rotor  data,  which  are  common  to  all  the  test  cases 
presented,  are  included  in  the  output  file  listing  (Appendix  J).  Values  of  the  new  input 
variables  for  the  new  models  are  included  in  the  'Main  Rotor'  subsection  of  the  'Input  Data' 
section  of  this  file. 
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Fig.  4  Effect  of  Wake  Model  Variations  on  Blade  Load  Distribution  for 
S-58  using  Vortex  Line  Model  (Ct/o  =  0.0817) 


Figs  4  and  5  show  the  effect  of  the  same  wake  model  variations  applied  to  the  vortex 
line  and  vortex  ring  models  respectively,  each  with  compressibility  included.  In  Table  1, 
the  maximum  blade  loading  (at  a  radial  station  of  0.95)  is  tabulated  for  these  variations,  but 
both  with  and  without  compressibility  included.  Each  model  gives  reasonably  similar 
distributions,  but  the  maximum  loading  given  by  the  vortex  line,  prescribed  wake  model 
with  a  concentrated  far  wake  is  about  10%  less  than  that  given  by  the  other  models. 


Fig.  5  Effect  of  Wake  Model  Variations  on  Blade  Load  Distribution  for 
S-58  using  Vortex  Ring  Model  (Ct/o  =  0.0317) 
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TABLE  1 

Effect  of  Wake  Model  Options  on  Maximum  Blade  Load  for  S-58 


Wake  Model 

Compressible 

Vortex  Line  Vortex  Ring 

Incompressible 

Vortex  Line  Vortex  Ring 

Free 

388 

375 

365 

353 

Prescribed  (concentrated  far  wake) 

349 

380 

316 

338 

Prescribed  (distributed  far  wake) 

388 

385 

346 

346 

5.  CONCLUDING  REMARKS 

By  incorporating  the  simplified  hover  wake  models  described  here  into  CAMRAD,  the 
models  themselves  have  been  enhanced  by  allowing  access  to  compressibility  effects 
included  in  the  two-dimensional  airfoil  tables  used  by  CAMRAD.  In  addition,  the  range  of 
rotor  wake  models  in  CAMRAD  has  been  broadened  and  now  includes  a  free  wake  hover 
model  (either  vortex  line  or  vortex  ring)  that  does  not  have  the  convergence  problems  of  the 
existing  free  wake  model  in  CAMRAD  when  used  for  hover. 

In  deciding  which  of  the  simplified  models  to  use,  consideration  should  be  given  to 
maintaining  a  balance  between  the  degrees  of  approximation  used  within  the  wake  model 
itself  Specifically,  when  representing  vortices  by  just  straight  lines  in  all  wake  regions,  the 
added  complexity  of  a  free  wake  solution  may  not  be  warranted.  The  more  physically 
accurate  representation  by  rings  would  therefore  seem  to  be  more  consistent  with  a  free 
wake. 
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APPENDIX  A 
Computational  Procedure 

Initial  Step  -  Free  Wake  only 

For  r  =  r*[ij  at  i  =  set 

w(r)  =  K(T) 
u(r)  =  0 

where 

r  =  radial  displacement  coordinate  along  the  blade 

fAli]  =  r  at  mid-points  of  aerodynamic  segments 

M  =  number  of  aerodynamic  segments 

u,  w  =  net  radial  and  axial  induced  velocity  components 

X,  =  vertical  component  of  induced  velocity,  given  by  uniform  inflow 
method  based  on  momentum  theory  (from  CAMRAD) 

Iteration  Loop 

STEP  1  .  Free  Wake  only 

Calculate  bound  circulation  distribution  r(r)  along  the  blade  for  r  =  r^li]  at  i  =  1  ....M: 
r(r)  =  y  r  c(r)  Cl  (r,  a,  rM,ip) 

where 

c  =  blade  chord 
Cl  =  blade  lift  coefficient 

The  latter  is  interpolated  from  2-D  airfoil  tables  as  a  function  of  angle  of  attack  a  and 
Mach  number  r  M.jp,  with  o  given  by 

o(r)  =  e(r)  -  «(r) 

where 

e(r)  =  blade  pitch  angle  (specified) 

(t)(r)  =  blade  inflow  angle  (  =  w(r)/r ) 
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STEP  2 

Representing  the  near  wake  by  semi-infinite  line  vortices  at  r  =  r^Efi]  for  i  = 
calculate,  using  the  Biot-Savart  law,  the  induced  velocity  w^fr)  along  the  blade  at  r  =  r^li] 
for  i  =  1,...M: 


Mtl 

w,(r)  =  -  ^  (r*Eli]  -  r  +  r^nli  +  r) 

i  •  1  ^  ^ 

where 

Ar(rAE[i|)  =  r(rA[i])  -  r(rA[i-l])  for  i  =  2,...M 

Ar(rAE[ll)  =  r(rA[ll) 

Ar(rAE[M-t-l])=  -  r(rA[M]) 

The  above  bound  circulation  distribution  is  given  in  Step  1  for  the  free  wake,  and  by  the 
uniform  inflow  method  based  on  momentum  theory  (from  CAMRAD)  for  the  prescribed 
wake. 


STEP  3 

Determine  the  circulation  strength  r,  and  location  r,  of  each  rolled-up,  concentrated 
vortex  along  the  blade  for  s  =  1„,.S  (tip  to  root),  where  S  is  the  total  number  of  these 
vortices.  The  boundaries  for  each  region  to  be  represented  by  concentrated  vortices  are 
first  defined  as  follows. 

The  outboard  boundary  of  the  tip  vortex  region  is  at  TAEfMol.  where  Mo  =  M  -i-i ,  and  the 
inboard  boundary  is  at  rAslM,],  where  .r(rA(M,])  is  the  maximum  of  the  circulation 
strengths  (r^,,)  calculated  in  Step  1  along  the  blade  at  the  aerodynamic  segment  mid¬ 
points. 

The  boundaries  for  the  inboard  regions  are  defined  to  be  at  the  closest  aerodynamic 
segment  boundary  inboard  of  the  values  defined  by  the  user.  The  array  indices  are  defined 
as  M2,...Ms.i,  corresponding  to  boundaries  rAE(M2],.„rAE[Ms.i].  The  most  inboard 
boundary  is  at  the  blade  root,  where  the  index  is  Ms,  and  the  boundary  r^ElMs]. 


a)  Free  Wake 

The  circulation  strength  r,  and  location  r,  (centroid  of  circulation  distribution  over  the 
region  r^ElMs]  torAEfM..,]  )  are  given  by 


r.  = 


X  Ar(rAE[i]) 

i.  M, 


r.  = 


r. 


rAE[il  Ar(rAE|i]) 


For  s  =  2,...S,  r,  is  then  scaled  by  an  empirical  factor  fr  (default  value  of  0.6). 
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b)  Prescribed  Wake 

Here  S  =  2,  and  the  concentrated  vortices  are  assumed  to  be  at  the  boundary  extremities 
(tip  and  root),  with  the  magnitude  of  the  circulation  strength  of  each  equal  to  r„,„,  i.e. 

n  =  I'm.!  r2  =  -  r„,„ 

ti  =  r*£[M+l]  rj  =  r,^[l] 

STEP  4  -  Free  Wake  only 

By  interpolating  induced  velocity  components  along  the  blade  calculated  at  the  previous 
iteration  in  Step  8  (zero  initially),  determine  values  at  each  concentrated  vortex  position,  for 
s  =  1,...S: 

{uil,  =  u,  at  r  =  r. 

=  u,(rA[i-l))  +  (u,(rA[i])-u,(rAli-l]))j^p^^|^jp^jj^j  for  rJi-1]  <  r.  < 

and  similarly  for  {w,}„  (uf).,  and  |wf),. 

STEP  S  •  Free  Wake  only 

Set  intermediate  wake  geometry,  defining  radial  and  axial  positions  of  the  concentrated 
vortices  at  each  vortex  line  or  ring  level,  i.e.  r,,  and  z,,„  for  s  =  1,...S  and  t  =  1,..,T, 
where  T  is  the  number  of  levels: 


r„  =  r„.,  +  AT.,, 


z,,  =  Z,.,.,  +  AZ,, 

where  the  incremental  displacements  Ar,.,  and  Az.,,  are  given  from  the  previous  iteration 
(Step  1 3),  but  are  approximated  initially  by 

Ar.,  =  0 

AZ.,  =  ^  X.(r*[M]) 

N 


STEP  6  -  Free  Wake  only 

Calculate  radial  and  axial  components  of  the  induced  velocity  at  each  vortex  position  in 
the  intermediate  wake,  i.e.  for  s  ^  1,...S  and  t  =  1,...T,  (a)  due  to  the  intermediate  wake 
to  give  (Ui)„  and  (w,).,,  and  (b)  due  to  the  far  wake  to  give  {uf}„  and  {wf1„  .  Only  the 
radial  components  are  shown  below;  the  axial  components  are  given  by  substituting  w  for  u 
in  all  expressions: 


11 


T 

^m.o*  ^m.n  "  ^m) 

ft  -I  (#  0 


{UfI 


a 

■  .t  ~  ^m.T+l  "  ^•.t»  ^m) 


Expressions  for  the  above  right-hand  side  velocity  components,  together  with 
equivalent  axial  components,  are  given  in  Appendix  B  for  both  vortex  line  and  vortex  ring 
models. 


STEP  7  -  Free  Wake  only 

Using  the  Biot-Savart  law  as  in  Step  2,  calculate  the  induced  velocity  (w,,),  at 
concentrated  vortices  on  the  blade  due  to  the  trailing  near  wake  for  s  =  1,...S  : 


{Wbl.  =  Wb(r)  atr  =  r. 


'  -1  <«•) 


STEP  8  -  Free  Wake  only 

Determine  net  radial  and  axial  components  of  the  induced  velocity  at  concentrated  vortex 
positions  on  the  blade  (i.e.  at  t  =  0 )  and  at  each  vortex  position  in  the  intermediate  wake; 


=  Iui  +  Uf), 

for 

t=0 

=  {Ui+Uf)„ 

for 

t  = 

1,...T 

=  {w,  +  w,.  + Wb).  +  {w«jf). 

for 

t=0 

=  (w,  +  Wpl.,,  +  {w,^|. 

for 

t  = 

1,...T 

The  self  induced  velocity  { w,^),  of  a  vonex  ring  of  radius  r,  is  given  by 


where 

rj  =  vortex  core  radius  r^  limited  to  a  minimum  of  0.005,  i.e.max(0.005,  rj 
ri  =  rix,  (burst  vortex  core  radius)  if  d^,  ^0  or  z, ,  <  d^, 

=  r„c  (unburst  vortex  core  radius)  otherwise 
dbv  =  core  burst  test  parameter  (<  0  to  suppress  bursting)  -  from  CAMRAD 
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STEP  9 
a)  Free  Wake 

Using  net  induced  velocity  components  at  the  blade  and  in  the  intermediate  wake, 
determine  new  wake  geometry  for  s  =  and  t  =  1,...T; 

'  r,^.i  +  Ar,i 

Z.,  =  +  AZ., 

where  the  incremental  displacements  ar,,  and  Az,,  are  given  by 


at,.,  = 


«(u.,,.i  -t-  u,.,) 
N 


AZ.,,  = 


ll(w,.,.i  +  w,,) 
N 


b)  Prescribed  Wake 

Using  prescribed  wake  geometry  based  on  either  of  the  options  in  CAMRAD,  i.e. 
(a)  Landgrebe  or  (b)  Kocurek  and  Tangier,  set  tip  and  root  vonex  positions  (noting  S  =  2) 
for  t  =  1,...T: 


r,,.  =  K.  +  (1  -K„)e-^’‘'^»'^ 
fit  =  r^Elll 

Z,,  =  Z2.  =  -  Ir 


where 

Ki,  Kj  =  axial  settling  rates  of  tip  vortex  before  and  after  passage  of  following  blade 
K3,  K4  =  radial  contraction  parameters  for  tip  vortex 

These  parameters  are  given  by  CAMRAD  on  specifying  the  appropriate  option,  i.e. 
value  of  OPRWG  in  Namelist  NLWAKE. 

STEP  10 

Calculate  radial  and  axial  components  of  the  induced  velocity  along  the  blade  (a)  due  to 
the  intermediate  wake  to  give  Ui(r)  and  w,(r),  and  (b)  due  to  t  ie  far  wake  to  give  upfr)  and 
WpCr)  for  r  =  rA[il  at  i  =  1,..,M.  As  in  Step  6,  only  the  radial  components  are  shown 
below,  with  axial  components  obtained  by  substituting  w  for  u; 

Ui(r)  =  ^  fd»  1  r7f] 

m  -I  „  ,1  V  '  y 
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for  distributed  far  wake 


s 

tipCf)  “  tip(r>  rm) 

m  >1 


-  Ui(r,  Tj,  Zi  7,  4ri)  +  U|(r,  r2,  Z2_T.  r'2) 


for  concentrated  far  wake 
(prescribed  only) 


where 


d'  =  (fm.^  -  r)^  +  (z„,  -  zy 

Expressions  for  the  above  right-hand  side  velocity  components,  together  with  the 
equivalent  axial  components,  given  in  Appendix  B  are  again  used. 

STEP  11 

Determine  net  radial  and  axial  components  of  the  induced  velocity  along  the  blade  to 
give  u(r)  and  w(r)  forr  =  rA[i]  at  i  = 

u(r)  =  u,(r)  -t-  UfCr) 
w(r)  =  w,(r)  -i-  wpfr)  Wb(r) 

For  prescribed  wake,  go  to  Step  14. 

STEP  12  -  Free  Wake  only 

Test  for  convergence  of  induced  velocity;  if  the  maximum  number  of  iterations  has  been 
reached  or 

M 

n  X  [  (w(rA[il))„..  -  {w(rAli]))„iJ  '  <  [j  w„„  e]  ' 

where 

=  max  lw(rA(i])l 

e  =  tolerance  for  induced  velocity  convergence 
and  then  go  to  Step  14. 

STEP  13  -  Free  Wake  only 

To  help  prevent  numerical  instability  in  the  iterative  procedure,  lag  new  (non- 
convergent)  solution  for  (a)  induced  velocity  u(r)  and  w(r)  along  the  blade  for  r  =  rA[i]  at  i 
=  1,...M,  and  (b)  wake  geometry  incremental  displacements  Ar,,  and  Az,,  for  s  =  1,.  .S 
andt=l,...T: 


u(r)  =  /  |u(r)|„,  +  (!-/)  (u(r)l,j„ 
w(r)  =  /  {w(r))„e.  +  (1  -/)  (w(r)),jj 
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Ar„  =  /  {Ar.,}^  +  (1  -/)  {Ar..,)„,d 
AZ.J.  =  /{^u)n~  +  (1  -/)  {AZ..,)old 

where  the  factor  /  used  to  introduce  lag  into  the  solution  is  typically  0.1. 

Having  completed  an  iterative  cycle  for  the  free  wake,  go  back  to  Step  1 . 

STEP  14 

Transform  induced  velocity  components  u(r)  and  w(r)  along  the  blade  to  longitudinal, 
lateral,  and  vertical  components  X,(r,  i((),  X,(r,  y),  and  X,(r,  y)  used  by  CAMRAD; 

X,  1  r-u  cosy 
X,  =  u  siny 
X,  w 
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APPENDIX  B 

Velocity  Components  Induced  by  Wake  Vortex  Elements 


Vortex  Line 


u,(r,  p,  h,  r)  -  ^  [  (p  +  ■  (p  -  r)2  +  h2  ] 

w,(r,  p,  h,  r)  =  ^  [  (p  +  r)2  +  h2  (p  -*r)2  +  1 

Wp(r.  p,  h.  r)  =  ^  [  11  -  arctan^^-^:]  * 

= 

=  [arctanj^^j  +  arctan^^^ 


Vortex  Ring 


u.(r,p,h,r)  = 

w.(r,p,h,r)=  [k - 

UF(r,  p.  h,  D  =  -  ^^lf  [  K(2  -  k^)  -  2E  ] 


E(2  -  k^)  1 

jrrwr  * 

E(1  -^kMl  +p/r))  j 


wp(r,  p,  h,  n  = 

n  '  ^ 


-  rcos^) 


2prcos<J> 


[ 


1  - 


for  p  <  r 
for  p  =  r 
for  p  >  r 


r^  +  p^  +  -  2prcos0 


d(|i 


where  the  latter  is  integrated  numerically,  and 


k2  = 


4pr 


(p  +  r)2  + 

E  =  1  +  i(F4)(l-k2)  +  ±(F-.||)(1  .k2)2 
K  =  F  +  KF-  1)(I  -k^)  +  ^(F-^)(I  -k^)2 
F  =  U^J==.] 
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APPENDIX  C 

Block  Diagram  for  Free  Wake  Method 


APPENDIX  D 

Block  Diagram  for  Prescribed  Wake  Method 
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APPENDIX  E 

Modified  CAMRAD  Subprograms 


SUBROUTINE  GE0KR1( LEVEL) 


CCOTMON  /KTIP/KT 


MOD 


CALCULATE  WAKE  GEOMETRY  DISTORTION 


FOR  KOCUREK  AND  TANGLER  AND  LANDGREBE  MODELS,  FWGT(l)  AND  FWGT(2)  MOD 
ARE  USED  AS  FACTORS  FOR  KT( 1)  AND  KT(2)  MOD 

KOCUREK  AND  TANGLER 
FB».  000729*™ 

FC-2.3-.206*™ 

FM-1.-.25*EXP(  .04*™) 

FN-.5-.0172*™ 

KT(l)-(FB+FC*(ABS(CTG)/FLQAT(NBLADE)**FN)**FM)*FWGTa)  MOD 

KT(2)-SQRT(ABS(CTG-FLOAT(NBLADE)**FN*(ABS(-FB/FC) )**(!. /FM) ) ) 
1*FWGT(2)  MOD 

KT(3)»4.*CTH 
KT(4)«.78 
GO  TO  17 
C  LANDGREBE 

16  KT(l)«.25*(CTOS+.001*TO)*FVgGT(l)  MOD 

KT(2)-(1.  +  .01*W)*CTH*FWGT{2)  MOD 

KT( 3)-. 145+27. *CTG 
KT(4)-.78 


SUBROUTINE  INPWl 

COMMON  /widata/fac™n,opvxvy,knw,krw,kfw,kj>i,rru,fru,pru,fnw,dvs,d 
1LS,C0RE(5) ,OPCORE(2) ,WKMODL(13) ,0PNWS(2),LHW,0PHW,0PRTS,VELB,DPHIB 
2 ,DBV,QDEBUG,MRG,NG( 30 ) ,MRL,NL( 30 ) ,OPWKBP( 3 ) , KRWG, OPRWG, FWGT( 2 ) , FWG 
351  ( 2 ) ,  FW3SO(  2 ) ,  KWGT(  4 ) ,  KWGSI  (  4 ) , KV4GSO(  4 ) 

INTEGER  OPVXVY , OPCORE , WKHODL , OPNWS , OPHW , OPRTS , OPWKBP , OPRWG 
1  ,NIVL,NIBV,WFMODL,OPMODL,ITERV  MOD 

REAL  KWGT,KWGSl,KWGSO 

1  ,RIBB(8) ,FGAHMA  MOD 

COMMON  /G1DATA/KFWG,OPFWG,ITERWG,FAC™g,WGMODL(2)  ,RTOG(2)  ,C0REWG(4 

1 ) , MRVBWG , LDMW3 , NDHWG ; 36 ) , I PWGDB ( 2 ) , QWGDB , DQWG ( 2 ) 

INTEGER  OPFW3,MGMODL 
COMMON  /TMnATA/IT«X(182) 

INTEGER  DEBUG 

EQUIVALENCE  (TMXX{ 41 ) .DEBUG) 

COMMON  /UNnNO/NFDAT,NFAFl,NFAF2,NFRS,NFEIG,NFSCR,NUDB,NUOUT,NUPP, 
INULIN.NUIN 

CaiMON  /RING/  NIBV,RIBB,NIVL,FACTlV,EPlVEL,WFMODL,OPMODL,FGAMMA,  MOD 
1  ITERV  MOD 
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READ  WAKE  NAMELIST 

NAMELIST  /NLWAKE/FACTWN,OPVXVY,KNW,KRW,KFW,KDW,RRU,FRU,PRU,FNW,DVS 

1 ,  DLS , CORE , OPCORE , WKMODL , OPNWS , LHW , OPHW , OPRTS , VELB , DPH I B , DBV , QDEBUG 

2 ,  MRG ,  NG ,  MRL ,  NL ,  OPWKBP ,  KRWG ,  OPRWG ,  FV«;T  ,  FWGSI ,  FWGSO ,  KWGT ,  KW3SI ,  KWGSO 

3 ,  KFWG , OPFWG , I TERWG , FACTWG , WGMODL , RTWG , COREWG , MRVBWG , LDHW3 , NDMM3 , 1 P 


4WGDB,C3WGDB,DQWG 

5  ,NIVL,RIBB,NIBV,FACTIV,EPIVEL,WFMODL,OPMODL,FGAMMA,ITEBV  MOD 

MOD 

- DEFAULT  VALUES  FOR  VORTEX  LINE/RING - MOD 

MOD 

NIBV=2  MOD 

DO  1=1,8  MOD 

RIBB(I)=0.0  MOD 

END  DO  MOD 

NIVL=4  MOD 

FACTIV=0 . 1  MOD 

EPIVEL=0.05  MOD 

WFMODL-2  MOD 

OPMODL=0  MOD 

FGAMMA=0 . 6  MC» 

ITERV=200  MOD 

- end - MOD 


999  FORMAT  (1X,33HREADING  NAMELIST  NLWAKE  (ROTOR  1)) 
WRITE  (NU0UT,999) 

READ  (NUIN, NLWAKE) 

IF  (DEBUG  .GE.  2)  WRITE  ( NUDB , NXWAKE ) 

RETURN 

END 


SUBROUTINE  PRNTWl 


COMMON  /R1DATA/R1XX(932)  MOD 

EQUIVALENCE  ( R1XX( 81 ) ,RROOT)  MOD 

EQUIVALENCE  (TMXX( 77 ) ,MPSI ) , (TMXX(157) ,LEVEL) 

COMMON  /UNITNO/NFDAT , NFAFl , NFAF2 , NFRS , NFEIG, NFSCR , NUDB , NUOUT , NUPP , 
INULIN.NUIN 

REAL  RIBB(8)  .FGAMMA  MOD 

INTEGER  NIVL,WFMODL,OPMODL,NIBV,ITERV  MOD 

COMMON  /RING/  NIBV,RIBB,NIVL, FACTIV, EPIVEL,WFHODL,OPMODL, FGAMMA,  MOD 
1  TTERV  MOD 

PRINT  WAKE  INPUT  DATA 


VORTEX  RING/LINE  -  - 

991  FORMAT (//IX, 'VORTEX  LINE  AND  VORTEX  RING  MODELS  (PRESCRIBED  AND  FR 

1  EE ) '/5X, 'NUMBER  OF  INTERMEDIATE  VORTEX  LEVELS,  NIVL  =' ,  I3/5X, 'FAR 

2  WAKE  MODEL  (0  TO  OMIT,  1  FOR  CONCENTRATED,  2  FOR  SHEET),  WFMODL 
3=' ,I3//5X, 'FOR  FREE  WAKE  ONLY ' //5X ,' FACTOR  INTRODUCING  LAG  IN  INDU 
4CED  VELOCITY,  FACTIV  =' ,F10.4/5X, 'TOLERANCE  FOR  INDUCED  VELOCITY, 
5  EPIVEL  =' ,F10.4/5X, 'ROLLED-UP  VORTEX  V/EIOTTING  FACTOR  (EXCLUDING 
6TIP),  FGAMMA  =' ,F9.4/5X, 'MAXIUM  NUMBL?.  OF  INDUCED  VELOCITY  ITERAT 
7IONS,  ITERV  =' ,I4/5X, 'NUMBER  OF  ROLLED-UP  VORTICIES,  NIBV  =',I3) 
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oo*‘Ooo  ....  ..  ••non 


986  FORMAT (/5X, 'INBOARD  EDGE  OF  ROLLED-UP  VORTICIES,  EXCLUDING  ROOT  AN 
ID  TIP  (NIBV-2  VALUES), '//15X,'RIBB  =',8(F9.4)) 

END - 


999  FORMAT  (//1X,23HN0NUNIF0RM  INFLOW  MODEL/  MOD 

*5X, 'VORTEX  LINE  MODEL  IF  1,  VORTEX  RING  MODEL  IF  2,  OPHODL  =',l3/MOD 
*  5X,27HEXTENT  OF  NEAR  WAKE, HOD 

1  KNW  =,I5/5X,33HEXTENT  OF  ROLLING  UP  WAKE,  KRW  = , I 5/5X , 26HEXTENT 

2  OF  FAR  WAKE,  KFW  = , I5/5X, 30HEXTENT  OF  DISTANT  WAKE,  KEW  =,I5/5X 
3,37HROLLUP  INITIAL  RADIAL  STATION,  RRU  =, FIO . 4/5X, 42HROLLUP  INITI 
4AL  TIP  VORTEX  FRACTION,  FRU  =, FIO. 4/5X, 27HROLLUP  EXTENT  (DEG),  P 
4RU  =,F10.2/5X,37HNEAR  WAKE  TIP  VORTEX  FRACTION,  FNW  = , FIO . 4/5X, 50 
5HNUMBER  OF  SPIRALS  IN  AXISYMMETRIC  FAR  WAKE,  LHW  =, I5/5X, 40HAXISY 
6MMETRIC  WAKE  GEOMETRY  IF  0,  OPHW  =,I3) 


WRITE  (NUOUT,990)  KRWG,OPRWG,  (  FWGT(  I ) ,  FWGSI(  I )  ,  FWGSO(  I )  ,  I=1,2),(K 


1WGT(I),KWGSI(I),KWGS0(I),  1=1,4) 

WRITE  (NUOUT,991)  NIVL.WFMODL,  FACTIV,  EPIVEL,  FGAMMA,  ITERV, NIBV  MOD 

IF  (NIBV  .GT.  2)  THEN  MOD 

IF  (RIBB(l)  .LT.  RROOT)  THEN  MOD 

DRI=(0.9-RROOT)/FLOAT(NIBV-1)  MOD 

DO  1=1, NIBV-2  MOD 

RIBB(I)=RROOT+I*DRI  MOD 

END  DO  MOD 

END  IF  MOD 

WRITE (NUOUT, 986)  (RIBB( I ) ,1=1 , NIBV-2)  MOD 

END  IF  MOD 

IF  (LEVEL  .LE.  1)  GO  TO  1 


SUBROUTINE  RAMF(LEVELl,LEVEL2,OPLMDA) 


COMMON  /RING/RIXX(16)  MOD 

INTEGER  OPRTR2,DEBUG,MHARH(2),MHARMF(2) 
l,OPMODL  MOD 

EQUIVALENCE  (TMXX( 77 )  ,MPSI ) ,  ( TMXX(  179 )  ,MHARM(  1 )  ) ,  (M 

1XX( 91 ) , ITERM ) , ( TMXX( 92 ) , EPMOTN ) , ( TMXX( 93 ) , ITERC ) , ( TMXX ( 94 ) , EPCIRC ) 

2 , ( TMXX( 49 ) , DEBUG) , ( THXX( 90 ) ,MREV) , ( TMXX( 89 ) ,MPSIR) , ( TRIMXX( 11 ) ,OPR 
3TR2  )  ,  (  BTR1XX(  4  )  ,CMEAN1 )  ,  (RTR1XX(  6 )  ,NBM1 )  ,  (  RTR1XX(  7  )  ,NTMl )  ,  (  RTRlXXC 
48 ) , NGMl ) , ( RTR2XX( 4 ) , CHEAN2 ) , ( RTF2XX( 6 ) , NBM2 ) , ( RrR2XX ( 7 ) , NTO2 ) , ( RTR 
52XX(8)  ,NGM2)  ,{B0DYXX(254),NAH),  (ENGNXX(  11 )  ,NDM)  ,  ( THXX(  181 )  ,MHARMF( 

61 ) ) , ( R1XX( 24 ) , SIGMAl ) , ( B2XX( 24 ) ,SIGMA2 ) 

7  , (RIXX(14) ,OPMODL)  MOD 

COMMON  /tlNlTNO/NFDAT, NFAFl , NFAF2 , NFRS , NFEIG ,  NFSCR , NUDB ,  NUOUT ,  NUPP , 
1NULIN,NUIN 

CALCULATE  ROTOR/AIRFRAME  PERIODIC  MOTION  AND  FORCES 


CALCULATE  INDUCED  VELOCITY 
CALL  WAKEUl 
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c - VORTEX  RINGAINE - 

IF(OPMODL.GT.O.AND.LEVELl.NE.O)  THEN 
CALL  PJAKERl  ( LEVELl ) 

GOTO  14 
ENDIF 

C - E34D - 

CALL  WAKENl  ( LEVELl ) 


MOD 

MOD 

MCO 

MOD 

M(X) 

MOD 

MOD 


C  END  MOTION  ITERATION 
C  TEST  CIRCULATION  CONVERGENCE 
OUT-0 

IF  (LEVELl  .EQ.  0)  GO  TO  53 

IF  ((LEVELl  .EQ.  2)  .AND.  (OPMODL  .GT.  0))  GO  TO  205  MOD 

G1MS=0. 


SUBROUTINE  TRIM 
COMMOJ  /TMDATA/rMXX(182) 

COMMON  /TRIMCM/TRIMXX(1604) 

COMMON  /CASECM/CASEXXO) 

COMMON  AING/RIXX(16)  MOD 

INTEGER  RESTRT,RSWRT,OPRTR2 

1  , OPMODL  MOD 

EQUIVALENCE  ( TMXX(  157),  LEVELl ) ,  ( ™XX(  158 ) ,  LEVEL2 ) ,  ( 'mXX(  159 ) ,  ITERU 

1)  , (TMXX(160) ,ITERR),(TMXX(161),ITERF),(TMXX(162),NPRNTT) , (TMXX(163 

2 )  , NPRNTP ) , ( TMXX ( 1 6 4 ) , NPRNTL ) , ( CASEXX ( 1) , RESTRT ) , ( CASEXX ( 5 ) , RSWRT ) , 
3(TRIMXX(11) ,OPRTR2) 

4  ,(RIXX(14), OPMODL)  MOD 

COMMON  AJNITNO/NFDAT , NFAFl , NFAF2 , NFRS , NFEIG , NFSCR , NUDB , NUOUT, NUPP , 
1NULIN,NUIN 

TRIM 


NONUNIFORM  INFLOW  AND  PRESCRIBED  WAKE 

2  IF  (LEVEL  .EQ.  1)  ITERR=MAX0( ITERR, 1 ) 

IF  (ITERR  .LE.  0)  GO  TO  3 

DO  20  IT=1, ITERR 

IF  (LEVELl  .EQ.  0)  GO  TO  22 

LEV1=1 

- VORTEX  LINE/RING  - 

- mod 

IF(OPMODL.GT.O)  THEN 

MOD 

CALL  GEOMRl(LEVl) 

MOD 

GOTO  22 

MOD 

END  IF 

MOD 

- END - 

- mod 

CALL  WAKECl(LEVl) 


C  NONUNIFORM  INFLOW  AND  FREE  WAKE 
3  ITERF»MAX0(ITERF,1) 

DO  30  IT-1, ITERF 
IF  (LEVELl  .EQ.  0)  GO  TO  32 
LEVl-LEVELl 
C 

IF(OPMODL.GT.O)  GO  TO  32  MOD 

C 

CALL  WAKECl(LEVl) 


SUBROUTINE  TIMER(N,I,T) 

COMMON  /TIMECM/rSTART( 14 ) ,TSUM( 14 ) ,NCALLS( 14 ) , IDBSAV( 23 ) , ICNT 
COMMON  /IMnATA/TMXX(182) 

EQUIVALENCE  (TMXX( 64 ) .DEBUG) , (TMXX( 41 ) , IDB( 1 ) ) , ( TMXX( 40 ) , ITDB ) 

INTEGER  DEBUG, IDB( 23) 

INTEGER*4  ITIME,HANDLE_ADR,CODE  MOD 

COMMON  /UNITNO/NFDAT ,  NFAFl ,  NFAF2 ,  NFRS ,  NFEIG ,  NFSCR ,  NUDB ,  NUOUT ,  NUPP , 
INULIN.NUIN 
C 

C  PROGRAM  TIMER 
C 

REAL  TFRACT(14) ,TCALL(14) 

INTEGER  ID(2,14) 

DATA  MT/14/ 

DATA  ID/4HCASE,4H  ,4HTRIM,4H  ,4HFLUT,4H  ,4HSTAB,4H  ,4H 
1TRAN,4H  ,4HSTAB,4HL  , 4HFLUT, 4HL  ,4HWAKE,4HC  ,4HGEOM,4HR 
2  ,4HRAMF,4H  ,4HMODE,4H  ,4HMOTN,4HR  ,4HPERF,4H  ,4HLOAD,4 
3H  / 

999  FORMAT  (1H1,17HC0MPUTATI0N  TIMES//50X,8HCPU  TIME, 4X,7HPERCENT,8X,6 
IHNUMBER , 3X , 9HTIME/CALL/52X, 5H( SEC ) , 19X , 8HOF  CALLS , 4X , 5H ( SEC ) / ) 

998  FORMAT  (1X,6HTIME  =,F12.3,4H  SEC) 

997  FORMAT  (1X,6HTIME  =,F12.3,4H  SEC, 5X,7H( START  ,2A4,1H)) 

996  FORMAT  (1X,6HTIME  =,F12.3,4H  SEC, 5X,5H( END  ,2A4,1H) ,5X,12H(CALL  NU 
1MBER,I3,18H,  TIME  INCREMENT  =,F12.3,5H  SEC)) 

901  FORMAT  (10X,4HCASE,31X,2F12.3,I12,F12.3) 

902  FORMAT  (lOX.llHTRlM  (TRIM) ,24X, 2F12 . 3 , 112 , F12 . 3 ) 

903  FORMAT  (10X,14HFLUTrER  ( FLUT) , 21X,2F12 . 3, 112, Fl2 . 3 ) 

904  FORMAT  (10X,22HFLIGHT  DYNAMICS  ( STAB) , 13X, 2F12 . 3 , 112, Fl2 . 3 ) 

905  FORMAT  ( lOX, 16HTRANSIENT  (TRAN) ,19X,2F12 . 3, 112, F12 . 3 ) 

906  FORMAT  ( 10X,23HLINEAR  ANALYSIS  (STABL) , 12X, 2F12 . 3 , I12,F12 . 3 ) 

907  FORMAT  ( lOX, 23HLINEAR  ANALYSIS  ( FLUTL) , 12X, 2F12 . 3 , 112 , F12 . 3 ) 

908  FORMAT  ( lOX, 25HNONUNIFORM  INFLOW  (WAKEC) , lOX, 2F12 . 3 , 112 , Fl2 . 3 ) 

909  FORMAT  (10X,21HWAKE  GEOMETRY  (GEOMR) ,14x, 2F12 . 3 , 112, Fl2 . 3 ) 

910  FORMAT  ( lOX, 25HVIBRATORY  SOLUTION  ( RAMF) , lOX, 2F12 . 3 , 112 , F12 . 3 ) 

911  FORMAT  (10X,18HROTOR  MODES  (MODE) ,17X, 2F12 . 3, 112 , F12 . 3 ) 

912  FORMAT  (10X,23HROTOR  EQUATIONS  (MOTOR) , 12X,2F12 . 3 , 112 , Fl2. 3 ) 

913  FORMAT  (10X,18HPERFORMANCE  ( PERF) ,  17X, 2F12 . 3 , 112 ,  n2 . 3 ) 

914  FORMAT  (10X,12HLOADS  (LOAD) ,23X,2F12. 3, 112, Fl2 . 3) 

CODE-2 

IF  (N  .EQ.  0)  GO  TO  10 
IF  (N  .EQ.  1)  GO  TO  11 
IF  (N  .EQ.  2)  GO  TO  12 
IF  (N  .EQ.  3)  GO  TO  13 
C  RETURN  PRESENT  TIME 

IERROR-LIB$STAT_TIMER  ( CODE , ITIME , HANDLEADR )  VAX 

T=.01*FLOAT(ITIME)  VAX 

IF  (DEBUG  .GE.  1)  WRITE  (NUOUT, 998)  T 
RETURN 

C  INITIALIZE 
10  CONTINUE 

IERROR=«LIB$INIT_TIMER  (HANDLE  ADR)  VAX 
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DO  100  JT=1,HT 
TSUM(JT)-0. 

NCALLS{JT)-0 
100  TSTART(JT)=0. 
rCNT=0 
RETURN 

C  START  TIMER 

11  CONTINUE 

IERROR*LIB$STAT_TIMER( CODE , ITIME, HANDLE_ADR )  VAX 

T-.01*FU)AT(ITIME)  VAX 

TSTART(I)=T 

IF  (DEBUG  .GE.  1)  WRITE  (NUOUT,997)  T, ID( 1 , 1 ) , ID( 2 , 1 ) 

IF  (I  .LE.  1)  GO  TO  113 
IF  (ICNT  .EQ.  1)  GO  TO  111 
tX3  112  11=1,23 
IDBSAV(II)=IDB(II) 

112  IDB(II)=0 
ICNT=1 

111  IF  (ITIME  .LT.  ITDB*1000)  GO  TO  113 
DO  114  11=1,23 
114  IDB(II)=IDBSAV(I1) 

113  CCWTINUE 
RETURN 

C  STOP  TIMER 

12  CajTINUE 

IERROR=LIB$STAT_TIMER(CODE, ITIME, HANDLE_ADR)  VAX 

T=.01*FLCIAT(  ITIME)  VAX 

DT=T-TSTART( I ) 

TSUM(I)-TSUM(I)+DT 
NCALLS ( I ) =NCALLS ( I ) +1 

IF  (DEBUG  .GE.  1)  WRITE  (NUOUT,996)  T,ID( 1, I ) , ID( 2, I) ,NCALLS( I ) ,DT 
RETURN 

C  PRINT  TIMES 

13  CONTINUE 
TCASE=TSUM(  1 ) 

IF  (TCASE  .NE.  0.)  TCASE=100 ./TCASE 

DO  130  JT=1,HT 

TFRACT ( JT ) =TSUM ( JT ) ‘TCASE 

TCALL(JT)=0. 

IF  (NCALLS(JT)  .NE.  0)  TCALL( JT)»TSUH(JT ) /FLOAT (NCALLS (JT) ) 

130  CONTINUE 

WRITE  (NUOUT,999) 

WRITE  (NUOUT,901)  TSUM( 1 ) ,TFRACT(1 ) ,NCALLS( 1 ) ,TCALL( 1 ) 

WRITE  (NUOUT,902)  TSUM( 2 ) ,TFRACT( 2 ) ,NCALLS( 2 ) ,TCALL( 2 ) 

WRITE  (NUOUT,903)  TSUM( 3) ,TFRACT( 3) ,NCALLS( 3) ,TCALL( 3 ) 

WRITE  (NUOUT,904)  TSUM( 4 ) ,TFRACT( 4 ) ,NCALLS( 4 ) , TCALL( 4 ) 

WRITE  (NUOUT,903)  TSUM(5),TFRACT(5),NCALLS(5),TCALL(5) 

WRITE  (NUOUT,906)  TSUM(6),TFRACT(6),NCALLS(6),TCALL(6) 

WRITE  (NUOUT,907)  TSUM(7) ,TFRACT(7 ) ,NCALLS( 7) ,TCALL( 7 ) 

WRITF  (NUOUT,908)  TSUM( 8 ) ,TFRACT( 8) ,NCALLS( 8 ) ,TCALL( 8 ) 

WRITE  (NUOUT,909)  TSUM(9 ) ,TFRACT(9) ,NCALLS( 9) ,TCALL( 9 ) 

WRITE  (NUOUT,910)  TSUM( 10 ) ,TFRACT( 10) ,NCALLS( 13 ) ,TCALL( 10 ) 

WRITE  (NUOUT,911)  TSUM( 11 ) ,TFRACT( 11 ) ,NCALLS( 11 ) , TCALL( 11 ) 

WRITE  (NUOUT,912)  TSUM( 12) ,TFRACT( 12) ,NCALLS( 12) ,TCALL( 12) 

WRITE  (NUOUT,913)  TSUM( 13 ) ,TFRACT( 13 ) ,NCALLS( 13 ) ,TCALL( 13 ) 

WRITE  (NUOUT,914)  TSUH( 14 ) ,TFRACr(14 ) ,NCALLS( 14 ) ,TCALL( 14 ) 

RETURN 

END 


n  n  n  n  no  nnnnnnnn 


APPENDIX  F 

Added  CAMRAD  Subprograms 


SUBROUTINE  WAKERK LEVEL) 


DETERMINE  THE  INDUCED  VEIOCITY  AT  THE  ROTOR  USING  A  VORTEX 
LINE  OR  RING  MODEL 

THE  WAKE  MAY  BE  EITHER  PRESCRIBED  OR  FREE 

THE  FAR  WAKE  MAY  BE  EITHER  NEGLECTED  OR  REPRESENTED  AS  A 
DISTRIBUTED  OR  CONCENTRATED  VORTEX 

INTEGER  DEBUGG,DEBUGV,I,IR,1T,J,JPSI,JR,JT,LEVEL,M,MPSI,MRA,N, 

1  NBLADE,NIBV,NTM,NTMl,OPCOMP,OPMODL,S,SMAX,T,NIVL,WFMODL, 

2  MS(0;10),ICOUNT,ITERV 

REAL  ALPHA ,  AREA ,  BETAC ,  BETAS ,  CL ,  COREB ,  CVERT ,  D5 ,  DA ,  DBV ,  DELW ,  DGAM , 

1  DGAMR ,  DU ,  CW ,  E ,  EP I VEL ,  EPR ,  ERR ,  FACTI V ,  FACTOR ,  FOLD ,  GAMR ,  H , 

2  LAMBDA, MACH, MTIP, PI, Ql,RMN,RROOT,T75,URES,WMAX,2ST,ANGL( 30) , 

3  DR(10,36) ,DZ( 10,36 ),GAMA( 30), GS( 10), R( 10,36 ) ,DROLD( 10,36 ) , 

4  RS(10),RUB(0:10),U(30),UIF(10,36),UNW(10,36),UOLD(30),US(10) , 

5  W(30),WB(30),WIF(10,36),WIFR(30) ,WNW( 10, 36 ) ,W0LD( 30 ) ,WRU( 10 ) , 

6  WS(  10 )  ,WSELF( 10), 2(10,36 ),DZOLD( 10,36 ) ,RIBB(8) ,FGAMMA 

CHARACTER  CHAR*  2 

-  CAMRAD  COMMON  BLOCK  - 

COMMON  /R1DATA/R1XX(932) 

COMMON  /RTR1CM/RTR1XX(1070) 

COMMON  /CONTCM/CONTXX(32) 

COMMON  /WKV1CM/WKV1XX(8165) 

COMMON  /TMDATA/TMXX{182) 

COMMON  /TRIMCMARIMXX(1604) 

COMMON  /W1DATA/W1XX(126) 

COMMON  /QR1CM/QR1XX(1139) 

COMMON  /MD1CH/MD1XX(6773) 

COMMON  /AEMNCM/AEMNXX( 78 ) 

COMMON  /RING/  NIBV,RIBB,NIVL, FACTIV, EPIVEL,WFMODL,OPMODL , FGAMMA, 

1  ITERV 

REAL  RA(  30 )  , 'IWIST(  30 )  ,CHORD(  30 ) , VIND(  3,30,36),  GAMOLD(  30,36), 

1  SINPSI ( 36 ) ,COSPSI ( 36 ) ,THET2L( 36 ) ,DRA( 30 ) , RAE{ 31 ) ,GAM( 30,36), 

2  ZETA(5,30),P1(5), CRCOLD( 36 ) , CRC( 36 ) , CORE( 2 ) 

EQUIVALENCE  ( RlXX(  150 )  ,HRA) ,  ( RTR1XX(  20 ) , RA(1 )  )  ,  ( CONTXX(  1 ) , T75 )  , 

1  (R1XX( 272 ) ,TOIST( 1 ) ) , (R1XX( 302) ,THETZL( 1 ) ) , (WlXX( 1 ) , FACTOR) , 

2  (R1XX(182),CH0RD(1)),(QR1XX(24),GAM(1,1) ),(R1XX(81) ,RROOT), 

3  (WKV1XX(1120),V1ND(1,1,1)),(R1XX(26),NBLADE),(TMXX(77),MPSI), 

4  (TRIMXX(57),SINPSI(1)),(R1XX(80),OPCOHP),(RTR1XX(7),NTM), 

5  (TRIHXX(21),COSPSI(l)),(WlXX(13),CORE(l) ),(W1XX(40),DBV), 

6  (RTR1XX( 2 ) ,MTIP) , (WKV1XX( 4360) , LAMBDA) , (QR1XX( 1104 ) ,CRC(  1 ) ) , 

7  (gRlXX( 22 ) , BETAC) , (QR1XX( 23 ), BETAS ) , ( RTR1XX( 50 ) ,URA( 1 ) ) , 

8  ( R1XX( 151 ) ,RAE{ 1 ) ) , (WKV1XX( 4 ) ,GAMOLD{ 1,1)), ( TMXX( 53 ) ,DEBUGV) , 

9  (TMXX(63),DEBUGG),(MD1XX(6148),ZETA(1,1)),(AEMNXX(31),P1(1)), 

A  (WKV1XX(1084),CRCOLD(1)) 

CCmON  /UNI  TNO/NFDAT ,  NFAFl ,  NFAF2 ,  NFRS ,  NFEIG ,  NFSCR ,  NUDB ,  NUOUT ,  NUPP , 
1NULIN,NUIN 


END  CAMRAD  COMMON  BLOCK 
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COMMOJ  /HELIC0M/Q1,W,U,GAMA,RS,GS,WS,US,DZ,R,DR,Z,WIFR,SMAX 


899 

900 

901 

902 

903 

904 


FORMAT (/IX, 'RING/LINE  LEVEL  ’,10(19)) 

FORMAT( /IX, 'VORTEX  NO. ',14,'  R=’ ,10(F9.5) ) 

FORMATdSX, 'Z=',10(F9.5))  inn 

F0RMAT(/1X, 'STRENGTH  OF  ROLLED  UP  VORTEX  ,F10. 


FORMAT(5X,10(F9.5) ) 
foRMAT(lX,A2,2X,10(F9.5) ) 


6) 


C 

C 

C 


INITIALIZE  VARIABLES 

data  PI/3. 14159265/ 

IF  (LEVEL  .EQ.  1)  THEN 
SMAX=2 
ELSE 

SMAX=NIBV 


END  IF 

MS(0)=HRA+1 

RUB(0)=RAE(MRA+1) 

MS(SMAX)=1 
RUB(SMAX)-=RB00T 
IF  (SMAX  .LE.  2)  GO  TO  20 
C  TEST  DATA  FOR  FATAL  ERRORS 
DO  I=1,NIBV-2 

IF  (I  .EQ.  1)  THEN 

IF  (RIBB(I)  .LT.  RROOT)  GO  TO  10 


10 


20 


ELSE 


TF  (RTBB(I)  .LT.  RIBB(I-l))  THEN 

WRITE(NUDB,*)  'ERROR  IN  DATA  : LOCATION  OF  INBOARD  VORTEX 


1  BOUNDARIES' 

GO  TO  10 
END  IF 
END  IF 
END  DO 

DO  I=l,NIBV-2 


S=SMAX-I 

RUB(S)=RIBB(I) 

END  DO 
GO  TO  20 

DBI= ( 0 . 9-RROOT )/FLOAT( NIBV-1 ) 
DO  I=0,NIBV-2 


S=SHAX-I 

RUB(S-1)=RUB(S)+DRI 


END  DO 

NTMl=MAX0(l,NTn) 


ICOUNT=l 

CVERT=180 ./PI 

Ql=2 . 0*PI/FLOAT( NBLADE ) 

IT=0 


C 

C  GEOMETRIC  PITCH  AND  INDUCED  VELOCITY  FROM  CAMRAD 
C  (FREE  WAKE  ONLY) 

C 

IF  (LEVEL  .EQ.  2)  THEN 
DO  IR=1,MRA 

ANGL( IR)=T75+(TWIST( IR)+THETZL( IR) )/CVERT 
W(IR)=-VIND(3,IR,1) 

DO  JT-1,NTM1 

ANGL ( IR ) =ANGL ( IR ) +ZETA ( JT, IR ) *P1 ( JT ) 

END  DO 
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END  DO 
END  IF 

POSITICWS  OF  ROLL-UP  BOUNDARIES 

IF  (SMAX  .GT.  2)  THEN 
DO  I=2,SMAX-1 
S=I 

DO  JR=1,MRA 
J=.JR 

IF  (RAE(J)  .GE.  (RUB(S)-l.E-4) )  GO  TO  101 

RUB(S)  IS  INBOARD  LIMIT  OF  THE  VORTEX  ROLL-UP  BOUNDARY 

END  DO 
)1  MS(S)=J 

END  DO 
END  IF 
MS(0)=MRA+1 
RUB(0)=RAE(MRA+1) 

MS(SMAX)-1 
RUB ( SMAX )=RROOT 

BEGINNING  OF  LOOP  FOR  NEXT  ITERATICN 

COMPUTE  BLADE  BOUND  CIRCULATIWJ 

D  IF( LEVEL  .EQ.  1)  THEN 

FROM  CAMRAD 

FOLD-1. O-FACTOR 
DO  JPSI=1,MPSI 

CRCOLD( JPSI )=FOLD*CRCOLD( JPSI )+FACTOR*CRC( JPSI ) 
DO  IR=1,MRA 

GAHOLD( IR, JPSI )=F0LD*GAM0LD( IR, JPSI )+FACTOR 
1  *GAM(IR,JPSI) 

GAMA( IR)-GAM0LD( IR, JPSI ) 

END  DO 
END  DO 
ELSE 

FROM  INDUCED  VELOCITY 
DO  80  JR-1,MRA 

DETERMINE  ALPHA, URES, MACH, COSL 

ALPHA- ( ANGL ( JR )-W( JR )/RA( JR ) ) *CVERT 
URES=W( JR ) *W( JR ) +RA( JR ) *RA( JR ) 

IF(URES.NE.O.O)  URES-SQRT(ABS(URES) ) 

MACH-URES*MTIP 

IF(OPCOMP.EQ.O)  HACH-0.0 

C  LIFT  COEFFICIENT 

CALL  AEROTl(ALPHA,MACH,RA(JR),l,CL,CD,CM) 

C  CIRCULATION 

GAMA( JR)=0 . 5*CL*URES*CHORD( JR) 
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80  CONTINUE 
END  IF 

LOOP  TO  FIND  STATION  OF  MAX.  CIRCULATION 

DO  90  JR=MRA,2,-1 
J=JR 

IF(GAHA(J)  .GT.  GAMA(J-l))  GO  TO  100 
90  CCOTINUE 

100  MS(1)-J  1  MS(1)  IS  LOCATION  OF  MAX.  CIRCULATICN  ON  BLADE 

RUB(1)=RAE(J) 

COMPUTE  INDUCED  VELOCITi  AT  BLADE  DUE  TO  NEAR  TRAILING  WAKE 
{DUE  TO  SHEET} 

105  DO  111  IR=1,MRA 
WB(IR)-0.0 
DO  110  JR=1,MRA+1 

IF  (JR  .EQ.  1)  THEN 

DELW=-l/( 4 . *PI ) * ( GAMA( JR ) ) * 

1  (1.0/(RAE(JR)-RA(IR))+1.0/(RAE(JR)+RA(IR) ) ) 

ELSE  IF  (JR  .EQ.  (MRA+1))  THEN 
DELW=1/(4.*PI)*(GAMA(JR-1) )• 

1  (1.0/(RAE(JR)-RA(IR))+1.0/(RAE(JR)+RA(IR) ) ) 

ELSE 

DELW=-l/( 4 . *PI ) * ( GAMA( JR ) -GAMA( JR-1 ) ) * 

1  (1.0/(RAE(JR)-RA(IR))+1.0/(RAE(JR)+RA(IR) ) ) 

END  IF 

WB(IR)=WB(IR)+DELW 

110  CONTINUE 

111  CONTINUE 

PRESCRIBED  WAKE  GEOMETRY 

IF( LEVEL  .EQ.  1)  THEN 
CALL  PRESWG(MS(1)) 

GO  TO  306 
END  IF 

FIND  CENTROID  OF  VORTEX,  RS(S) 


DO  S=1,SMAX 
GS(S)=0.0 
GAMR=0.0 

DO  120  JR=MS(S)+1,MS(S-1) 

IF  (JR  .EQ.  2)  THEN 
DGAM=GAMA(JR) 

ELSE  IF  (JR  .EQ.  (MRA+1))  THEN 
DGAH=-GAMA(JR-1) 

ELSE 

DGAM=GAMA( JR)-GAMA( JR-1 ) 

END  IF 

GS(S)=GS(S)-DGAM 

GS(S)  IS  THE  CIRCULATION  OF  THE  VORTEX,  S 

DGAMR=DGAM*RAE(  JR) 

GAMR=>GAMR+DGAMR 
120  CONTINUE 

IF  (ABS(GS(S))  .LT.  l.E-10)  THEN 
RS( S ) =( RUB( S )+RUB( S-1 ) )/2 . 0 
ELSE 
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RS(S)=ABS(GAMR/GS(S))  !  LOCATION  OF  VORTEX,  S 
IF  (S  .NE.  1)  GS(S)=GS(S)*FGAHMA 
END  IF 
END  DO 
C 

C  FIND  VELOCITIES  AT  RS(S),  WS(S),  AND  US(S) 

C 

DO  I=1,SMAX 
S=I 

DO  J=MS(S),MS(S-1) 

JR=J 

IF  (RS(S)  .LT.  RAE(aR))  GO  TO  160 
END  DO 
C 

C  INDUCED  VELOCITY 
C 

160  CALL  IVTERP(S,JR) 

C 

IF  (IT  .EQ.  0)  THEN 
WS(S)=0.0 
US(S)=0.0 
END  IF 
END  DO 

IF  (IT  .EQ.  2)  THEN 
DO  S=1,SMAX 
DO  T=1,TMAX 

IF  (T  .EQ.  1)  THEN 
Z(S,T)=DZ(S,T) 

R(S,T)=RS(S)+DR(S,T) 

ELSE 

Z(S,T)»Z(S,T-1)+DZ(S,T) 

R(S,T)=R(S,T-1)+DR(S,T) 

END  IF 
END  DO 
END  DO 
ELSE 
C 

C  ESTABLISH  INITIAL  WAKE  GEOMETRY  FROM  MOMENTUM  THEORY 
C 

DO  221  S=1,SMAX 
DO  220  T=1,NIVL 
IR=MS(1) 

219  D5=W(IR)*Q1 

C  CHECK  TO  SEE  IF  INFLOW  AT  STATION  IR  IS  .NE.  TO  ZERO. 

C  IF  INFLOW  AT  IR  .EQ.  ZERO  THEN  MOVE  IN  TO  NEXT  RADIAL  STATION. 

C  THIS  PREVENTS  PROGRAM  FROM  CRASHING  S,  HELPS  CONVERGENCE. 

IF  (D5  .EQ.  0.0)  THEN 
IR=IR+1 
GO  TO  219 
ENDIF 

DZ(S,T)=D5 

R(S,T)=RS(S)  !  INITIALISE  RADIAL  VORTEX  POSITONS 

DR(S,T)=0.0 

IF  (T  .EQ.  1)  THEN 

Z{S,T)=DZ(S,T)  !  INITIALISE  VERTICAL  VORTEX  POSITONS 
ELSE 

Z(S,T)=Z(S,T-1)+DZ(S,T) 

END  IF 

220  CONTINUE 

221  CONTINUE 
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END  IF 
C 

C  COMPUTE  VELOCITIES  IN  WAKE  DUE  TO  INTERMEDIATE  AND  FAR  WAKES 
C 

230  DO  271  S=1,SMAX 

DO  270  T=1,NIVL 

DZOLD(S,T)=D2(S,T) 

DROLD(S,T)=DR(S,T) 

C 

C  CALCULATE  VELOCITES  AT  R(S,T).  Z(S,T)  DUE  TO  INTERMEDIATE  AND  FAR  WAKES 
C 

CALL  VTXIF(R(S,T),Z(S,T),0,W1F(S,T),UIF(S,T) ) 

C 

C 

C  VELOCITY  IN  WAKE  DUE  TO  ROLLED-UP  NEAR  WAKE 
C 

H=-Z(S,T) 

IF  (OPMODL  .EQ.  1)  GO  TO  260 
WNW(S,T)=0.0 
UNW(S,T)=0.0 
DO  M=1,SHAX 
RHN=RS(M) 

CALL  IRING(R(S,T),RMN,H,GS(M),DU,EW) 

WNW(  S , T )  =WNW(  S  ,T ) +DW 
UNW(S,T)=UNW(S,T)+DU 
END  DO 
GO  TO  666 

260  WNW(S,T)=0.0 

UNW(S,T)=0.0 
DO  H=1,SMAX 

CALL  ILINE(R(S,T) ,RS(M) ,H,GS(M) ,DU,CW) 

WNW(S,T)=WNW(S,T)+DW 
UNW( S , T ) =UNW( S , T) +DU 
END  DO 

666  COREB*CORE(l) 

IF  (DBV  .GE.  0.0  .AND.  Z(l,l)  .LT.  DBV)  COREB=CORE ( 2 ) 

IF  (COREB  .LT.  0.005)  COREB=0.005 
C 

C  SELF  INDUCED  VELOCITY 
C 

WSELF(S)=1/(4.*PI)*GS(S)*(LCIG(8.0*RS(S)/COREB)-0.25) 

WIF(S,T)=WIF(S,T)+WNW(S,T)*0.5+WSELF(S) 

UIF{S,T)=UIF(S,T)+UNW(S,T)*0.5 

270  CONTINUE 

271  CONTINUE 
C 

C  COMPUTE  VELOCITY  AT  A  ROLLED-UP  NEAR  WAKE  DUE  TO  OTHER  ROLLED-UP 
C  NEAR  WAKES 
C 

DO  305  S=1,SHAX 
WRU(S)=0.0 
DO  M=1,SMAX 

IF  (RS(S)  .GT.  (RS(H)-l.E-4)  .AND.  RS(S)  .LT.  (RS(M) +1 . E-4 ) 

1)  THEN 

DW=GS(M)/(2.0*RS(M) ) 

ELSE 

EW=GS(M)*(1.0/(RS(M)-RS(S))+1.0/(RS(H)+RS(S))) 

END  IF 

WRU(S)=WRU(S)+CW 
END  DO 
C 

C  COMPUTE  NEW  WAKE  GEOMETRY 
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c 

Z(S,l)-( (WS(S)+WIF{S,1)+WRU(S)*1/(4.*PI))*0.5+WSELF(S)*0.5)*Q1 

DR(S,1)-(US(S)+UIF(S,1))*0.5*Q1 

DZ(S,1)-Z(S,1) 

R(S,1)=DR(S,1)+RS(S) 

290  DO  300  T=2,NIVL 

DZ(S,T)=Q1*(WIF(S,T)+WIF(S,T-1))*0.5 

Z(S,T)=Z(S,T-1)+DZ(S,T) 

DR(S,T)-Q1*(UIF(S,T)+UIF(S,T-1) )*0.5 
R(S,T)-R(S,T-1)+DR(S,T) 

300  CONTINUE 

305  CONTINUE 

Q  *  A  A  A  4r  *  -Ar  A  *  *  *  *  *  A  A  *  *  4^  *  *  *  A  *  *  *  *  *  4r  *  *  *  *  *  4r  A  *  *  *  *  *  *  ^  it 

IF  (DEBUGG  .GE.  2)  THEN 
WRITE(NUDB,*)  '  ' 

WRITE(NUDB,*)'  VORTEX  LINE/RING  GEOMETRY' 

WRITE(NUDB,*) '  ITERATICW  NUMBER' , ICOUNT 
WRITE (NUDB, 899)  (T,T=0,NIVL) 

DO  S«1,SMAX 
ZS=0 . 0 

WRITE (NUDB, 900)  S,RS(S) , (R(S,T) ,T=1 ,NIVL) 

WRITE(NUDB,901)  ZS, (Z( S,T) ,T=1 ,NIVL) 

WRITE( NUDB, 902)  GS(S) 

WRI'rE(NUDB,*)  '  ' 

END  DO 
END  IF 

Q**A**A**A* AAAAA AAAAAAAAAAAAAAAAAAAAA******************* 

C 

C  COMPUTE  NEW  INDUCED  VELOCITIES  AT  ROTOR 
C 

WMAX-0.0 

ERR-0.0 

306  DO  330  IR-1,MRA 
ZST-0.0 

UOLD(IR)-U{IR) 

WOLD(IR)=W(IR) 

IF  (IT  .EQ.  0)  UOLD(IR)=0.0 
C 

C  VELOCITIES  AT  RA{IR)  DUE  TO  INTERMEDIATE  AND  FAR  WAKES 
C 

CALL  VTXIF(RA(IR),ZST, LEVEL, WIFR(IR),U(IR)  ) 

C 

W( IR)=WIFR( IR)+WB( IR) 

IF  ((W(IR)*W(IR))  .GT.  WMAX)  WMAX=(W(  IR)  *W(  IR)  ) 

ERR-ERR+(W( IR)-M0LD( IR) )**2 
330  CONTINUE 
C 

C  NO  LOOP  IF  PESCRIBED  WAKE 
C 

IF  (LEVEL  .EQ.  1)  THEN 
IT-3 

GO  TO  380 
END  IF 
C 

C  COUNT  ITERATIONS  AND  CHECK  CONVERGENCE 
C 

ICOUNT-ICOUNT+1 
ERR— ERR/MRA/MRA 
EPR-EPIVEL*EPIVEL*WHAX/4 . 0 
IF  (ERR  .GT.  EPR)  GO  TO  350 
IT-3 

GO  TO  380 
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WEIGHT  NEW  VELOCITIES  AND  DISPLACEMENTS  FX)R  NEXT  ITERATION 

350  DO  360  IR»1,MRA 

FOLD=l . 0-FACTIV 

W( IR)=FOLD*WOLD( IR)+FACTIV*W( IR) 

U( IR)=FOLD*UOLD( IR)+FACTIV*U( IR) 

360  CONTINUE 
IT=2 

DO  375  S=1,SMAX 
DO  370  T=0,NIVL 

DZ ( S , T ) =FOLD*DZOLD ( S , T ) +FACTIV*DZ ( S , T ) 

DR ( S , T ) -FOLD*DROLD( S , T ) +FACTIV*DR(  S , T ) 

370  CONTINUE 
375  CONTINUE 

IT=0  : FIRST  ITERATION,  IT=2  :PROCEEDING  ITERATIONS,  IT-3  :CONVERGED 

380  IF  (IT  .EQ.  3)  GO  TO  395 

IF  (ICOUNT  .GT.  ITERV)  GO  TO  396 
GO  TO  70 
C 

395  COTTINUE 

IF  (DEBUGV  .GE.  1)  THEN 

WRITE(NUDB,*)  '  INDUCED  VELOCITIES  FOR  VORTEX  LINE/RING  MODEL' 
WRITE(NUDB,*)  '  RADIAL  STATIONS' 

WRITE(NUDB,903)  (RA(IR),IR=1,MRA) 

WRITE(NUDB,*)  '  BLADE  AXES' 

WRITE(NUDB,*)  '  AXIAL  VELOCITY' 

WRITE(NUDB,903)  (W(IR),IR-1,MRA) 

WRITE(NUDB,*)  '  RADIAL  VELOCITY' 

WRITE(NUDB,903)  (U(IR),IR=1,MRA) 

WRITE(NUDB,*)  '  ' 

END  IF 

DO  JPSI=1,MPSI 
DO  1R-1,MRA 

VIND(  3 ,  IR,  JPSI )— W(  IR) 

VIND( 2 , IR , JPSI ) -U( IR ) *SINPSI ( JPSI ) 

VIND ( 1 , IR , JPSI ) =U( IR ) *COSPSI ( JPSI ) 

END  DO 

IF  (DEBUGV  .GE.  2)  THEN 

WRITE(NUDB,*)  '  SHAFT  AXES' 

PSI=FLOAT( JPSI ) *360 . 0/FLQAT( MPSI ) 

WRITE(NUDB,*)  '  PSI  =',PSI 
CHAR='LX' 

WRITE(NUDB,90<1)  char,  (VIND(1,IR,  JPSI  ),IR=1,MRA) 

CHAR='LY' 

WRITE(NUDB,90‘3)  CHAR,(VIND(2,IR,JPSI),IR=1,MRA) 

CHAR= ' LZ ' 

WRITE(NUDB,904)  CHAR, ( VIND( 3 , IR, JPSI ) , IR-1 ,MRA) 

WRITE(NUDB,*)  '  ' 

END  IF 
END  DO 

IF  (DEBUGG  .EQ.  1)  THEN 
WRITE(NUDB,*)  '  ' 

WRITE(NUDB,*) '  VORTEX  LINE/RING  WAKE  GEOMETRY' 

WRITE ( NUDB , * ) '  NUMBER  OF  ITERATIONS ' , ICOUNT 
WRITE(NUDB,899)  (T,T=0,NIVL) 

EX)  S=1,SHAX 
ZS-O.O 

WRITE (NUDB, 900)  S,RS(S) , (R(S,T) ,T=1,NIVL) 
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WRITE(NUDB,901) 
WRITE(NUDB,902) 
WRITE(NUDB,*)  ' 
END  DO 
END  IF 


ZS.(Z(S,T),T=1,NIVL) 

GS(S) 

t 


C  CALCULATE  MEAN  INDUCED  VELOCITY 
LAMBDA-0.0 
AREA-0.0 
DO  IR=1,MRA 

nA=RA(IR)*DRA(lR) 

AREA-AREA+DA 
DO  JPSI=1,MPSI 

LAMBDA=LAMBDA-(VIND(  3,IR,  JPSI)-BETAC*VIND(  1,  IR,  JPSI ) 

1  -BETAS*VIND(2,IR,JPSI))*nA 

END  DO 
END  DO 

LAMBDA=LAMBDA/( FLQAT( MPSI ) ‘AREA ) 

RETURN 

396  WRITE(NUDB,*)  SOLUTION  NOT  CONVERGING  (INDUCED  VELOCIT 

lY)  ******' 

GO  TO  395 


C 


END 


SUBROUTINE  VTXIF(RST,ZST, LEVEL, WIF.UIF) 


SUBROUTINE  TO  CALCULATE  THE  VELOCITY  AT  (RST,ZST)  DUE  TO 
THE  INTERMEDIATE  AND  FAR  WAKES 


REAL  UIF,WIF,UI,WI,RMN,H,i:W,DU,CORE(2),KT, 

1  COREB,DBV,RST,DS,FACT,WP,UF,RIBB(8) 

INTEGER  M , SMAX , N , NI VL , OPMODL , Wi  MODL , LEVEL 

COMMON  /W1DATA/W1XX(126) 

EQUIVALENCE  {WlXX( 13 ) ,CORE( 1 ) ) , (W1XX( 40 ) , DBV) 

COMMON  /RING/  NIBV, RIBB,NlVL,FACTIV,EPIVEL,WFMODL, OPMODL, FGAMMA, 
1  ITERV 

COMMON  /KTIP/KT(4) 

REAL  Ql,W(30),GAMA(30),RS(10),GS(10),WS(10),US(10),D2tlO,36), 

1  R(10,36),DR(10,36),Z(10,36),WIFR(30),U!30) 

COMMON  /HELIC0M/Q1,W,0,GAMA,RS,GS,WS,US,DZ,R,DR,Z,WIFR,SMAX 


UIF-0.0 
WIF-0.0 
DO  M=1,SMAX 
WI=0.0 
UI=0.0 


INTERMEDIATE  WAKE 

DO  N=1,N1VL 
RMN=R(M,N) 

H=2(M,N)-ZST 
IF  (OPMODL  .EQ.  1)  THEN 

CALL  ILINE(RST,RMN,H,GS(M) ,DU,rW)  !  VORTEX  LINE 
ELSE 
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c 

c 

c 


c 

c 

c 


c- 


CALL  IRING(RST,RMN,H,GS(M)  ,IXJ,CW)  !  VORTEX  RING 
END  IF 

EFFECT  OF  VISCOUS  CORE  USING  APPROXIMATION  BY  SCULLY 

IF  (LEVEL  .GE.  1)  THEN 
COREB=CORE( 1 ) 

IF  (DBV  .GT.  0.0  .AND.  DBV  .GT.  Z(l,l))  COREB=COBE ( 2 ) 
DS= ( RMN-RST ) *  *2+H*H 
FACT=DS/( DS+COREB*COREB ) 

ELSE 
FACT=1 . 0 
END  IF 

WI=WI+DW*FACT 
UI=UI+DU*FACT 
END  DO 

FAR  WAKE 

H=Z  ( H ,  NIVL )  -2ST+DZ  ( M ,  NIVL ) 

RMN=R  ( M ,  NI VL ) +DR  ( M ,  NI VL ) 

IF  (OPMODL  .EQ.  1)  THEN 

CALL  FLINE(RST,RMN,H,GS(M)  ,UF,WF)  !  VORTEX  LINE 
ELSE 

CALL  FRING(RST,RMN,H,GS(M),UF,WF)  !  VORTEX  RING 
END  IF 

IF  (LEVEL  .EQ.  1)  THEN  !  PRESCRIBED  WAKE 

IF  (WFMODL  .EQ.  1)  THEN  !  CONCENTRATED  FAR  WAKE 

IF  (M  .EQ.  1)  THEN  !  TIP  VORTEX 

GAMM=4.0*GS(M) 

IF  (OPMODL  .EQ.  1)  THEN 

CALL  ILINE(RST,RMN,H,GAHM,UF,WF)  !  VORTEX  t tnE 
ELSE 

CALL  IRING(RST,BMI’,H,GAHH,UF,WF>  !  VORTEX  RING 
END  IF 

ELSE  !  RC)OT  VORTEX 

IF  (0"M0DL  .EQ.  1)  THEN 

CALL  ILINE(RST,R«N,H,GS(M),UF,WF)  !  VORTEX  LINE 
ELSE 

CALL  IRING(RST,RMN,H,GS(M) ,UF,WF)  !  VORTEX  RING 
END  IF 
END  IF 
END  IF 

IF  (WFMODL  .EQ.  0)  THEN  !  NO  FAR  WAKE 

WIF=WIF+WI 
UIF=UIF+UI 

ELSE  IF  (WFMODL  .EQ.  1)  THEN  !  CONCENTRATED  FAR  WAKE 
WIF=WIF+(WI+WF) 

UIF=UIF+(UI+UF) 

ELSE 

WIF=WIF+(WI+WF/DZ(M,NIVL)  )  !  DISTRIBUTED  FAR  WAKE 

UIF=UIF+(UI+UFA>Z(M,NIVL)  ) 

END  IF 
ELSE 

IF  (ABS(DZ(M,NIVL) )  .LT.  0.0001)  DZ( M,NIVL )=0 . 1 
WI F=WI F+ ( WI +WF/DZ ( M , NI VL ) )  !  FREE  WAKE 

UIF=UIF+(UI+UF/DZ(H,NIVL) ) 

END  IF 
END  DO 
RETURN 
END 


34 


oooonno  nonnn  nn  ooon 


SUBROUTINE  IVTERP(S,I) 


SUBROUTINE  TO  CALCULATE  THE  INDUCED  VELOCITY  AT  RS(S)  BY 
INTERPOLATING  THE  VELOCITIES  AT  RA(I)  (I=1,MRA) 

REAL  RA(30),GRAD 

INTEGER  S,I 

COMMON  /RTR1CM/RTR1XX(1070) 

EQUIVALENCE  (RTR1XX(20)  ,RA(1 )  ) 

REAL  Q1,W(30) ,GAMA(30) ,RS(10),GS{10)  ,WS(10) ,US(10) ,DZ{10,36) , 
1  R(10,36),DR{10,36),Z(10,36),WIFR(30),U(30) 

COMMON  /HELICOM/Ql,W,U,GAMA,RS,GS,WS,US,DZ,R,DR,Z,WIFR,SMAX 


GRAD=(RS(S)-RA(I-1) )/(RA(I)-RA(I-l) ) 
VERTICAL  VELOCITY 

WS( S)=WIFR( I-l )+(WIFR( I )-WlFR( I-l ) ) *GRAD 
RADIAL  VELOCITY 

US(S)*U(I-1)+(U(I)-U(I-1) )*GRAD 

RETURN 

END 


SUBROUTINE  ILINE( R, P,H, GAMMA, UILINE.WILINE) 

FUNCTION  FOR  RADIAL  AND  AXIAL  COMPONENT  OF  VELOCITY  INDUCED 
BY  VORTEX  LINE 

REAL  R,P,H, GAMMA, UILINE,WILINE 
DATA  PI/3.14159265/ 

IF  (ABS(H)  .LT.  1.E4  .AND.  ABS(P-R)  .LT.  l.E-4)  THEN 
UILINE=0.0 
WILINE^O.O 
RETURN 

END  IF 

W1LINE=GAMHA/(2.*PI)*( (P-F)/((P-R)**2+H*H)+(P+R)/( (P+R)**2+H*H) ) 
UILINE=-GAMMA/(2.*PI)*(H/( (P'R)**2+H*H)-H/i (P+R)**2+H*H) ) 

RETURN 

END 


SUBROUTINE  IRING(R,P,H, GAMMA, UIRING, WIRING) 


FUNCTION  FOR  RADIAL  AND  AXIAL  COMPONENT  OF  VELOCITY 
INDUCED  BY  VORTEX  RING 

REAL  R,P,H,K2,E,K,TEMP2 

DATA  PI/3. 14159265/ 

IF  (ABS(H)  .LT.  l.E-4  .AND.  ABS(P-R)  .LT.  l.E-4)  THEN 
WIRING=0 . 0 
UIRING=0 . 0 
RETURN 
END  IF 

CALL  ELLIPC0N(R,P,H,K2,E,K) 
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C 

TEMP2=K2/(R*P) 

IF  (TEHP2  .LE.  0.0  )  THEN 

WRITE(6,*)  '  BAD  SI®J.  MODULE  IRING’ 

STOP 
END  IF 

C  AXIAL  VELOCITY 

WIRING-SQRT( TEHP2 ) * ( K-E* ( 1 . 0- . 5*K2* ( 1 . 0+P/R) )/( 1 . 0-K2 ) ) 
WIRING=GAMMA/( PI*4 . ) *WIRING 
C  RADIAL  VELOCITY 

UIRING-H/( 2 . 0*R ) *SQRT( TEMP2 ) * ( E* ( 2 . 0-K2 )/( 1 . 0-K2 ) -2 . 0*K ) 
UIRING — GAMMA/( 4 . *PI ) *UIRING 
C 

RETURN 

END 

C - 

SUBROUTINE  ELLIPCON(R,P,H,K2,E,K) 

EVALUATE  CCWSTANTS  USED  IN  INDUCED  VELOCITY  COMPCWENTS 
DEFINED  BY  ELLIPTIC  INTEGRALS 

REAL  R,P,H,K2,E,K,TEMP2,F 

K2«4.0*R*P/( (R+P)**2+H»H) 

IF  (K2  .GE.  (1.0-2.E-8))  GO  TO  100 

TEMP2=1 . 0-K2 
F=LOG( 4 . 0/SQRT( TEMP2 ) ) 

E=1.0+.5*(F-.5)*(1.0-K2)+3./16.*(F-13./12.)*(1.0-K2)**2 
K-F+.25*{F-1.0)*(1.0-K2)+9./64.*(F-7./6.)*(1.0-K2)**2 

RETURN 

100  E-1.0 
K=10.0 

K2-1.0-3.0E-8 

RETURN 
END 


SUBROUTINE  FRING( R, P,H,GAMMA,UFRING,WFRING) 

RADIAL  AND  AXIAL  COMPONENT  OF  VELOCITY  INDUCED  BY  SEMI¬ 
INFINITE  VORTEX  CYLINDER 

REAL  R,P,H,PI,P3,P4,X3,I4,K2,E,K,TEMP2 

DATA  PI/3.14159265/ 

IF(ABS(P)  .LT.  l.OE-6)  P=1.0E-6 

AXIAL  VELOCITY 

WFRING=0 . 0 
P3=0.2*PI 

DO  20  P4=P3*0.5,PI-P3*0.5,P3 
X3=P*P+R*R+H*H-2 . 0*R*P*COS( P4 ) 

I4=p3*P*(P-R*C0S(P4) )/(P*P+R*R-2.0*R*P*COS(P4) ) 

*  *(1.0-H/SQRT(X3) ) 

C 


WFRING=WFRING+I 4 
20  CONTINUE 
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WFFING=GAMMA/( 2 . *PI ) *WFRING 


RADIAL  VELOCITY 

CALL  ELLIPC0N(R,P,H,K2,E,K) 

TEMP2=P/RA2 

IF  (TEMP2  .LE.  0.0)  THEN 
WRITE(6,*)  P,R,K2 

WRITE(6,*)  'BAD  SIGN.  MODULE  FRING' 

STOP 
END  IF 

UFRING=-GAMMA/( 2 . *PI ) *SQRT(TEMP2 ) * ( K* ( 2 . 0-K2 ) -2 . 0*E ) 

RETURN 

END 


REAL  FUNCTION  FLINE{R,P,H, GAMMA, UFLINE,WFLINE) 

RADIAL  AND  AXIAL  VELOCITY  DUE  TO  VORTEX  SHEET 

REAL  R,P,H,PI,E,GAMMA,UFLINE,WFLINE 

DATA  PI/3.1415927/ 

AXIAL  VELOCITY 

WFLINE=PI/2.-ATAN{H/(R+P) ) 

IF{P  .GT.  {R+I.OE-8))  WFLINE=PI-ATAN(H/(P-R) )-ATAN(H/(P+R) ) 
IF{P  .LT.  (R-l.OE-8))  WFLINE=ATAN(H/(R-P) )-ATAN(H/(R+P) ) 
WFLINE=GAMMA/(  2  .  *PI )  *WFLINE 
RADIAL  VELOCITY 

UFLINE=-GAMMA/(  4  .  *PI)  *LOG(  (H»H+(R+P)  **2  )/(  H*H+(  P-R )  **2  )  ) 

RETURN 

END 


SUBROUTINE  PRESWG(I) 

DETERMINES  PESCRIBED  WAKE  GEOMETRY  FOR  VORriiX  RING  AND  LINE  METHODS 


REAL  HV,RT,RROOT,KT(4),RIBB(8) 

INTEGER  IR,NIVL,WFMODL 

C  -  CAMRAD  COMMON  BLOCKS  - 

COMMON  /R1DATA/R1XX(932) 

COMMON  /KTIP/KT 
EQUIVALENCE  ( RlXX( 81 ) ,RROOT) 

C - END - 

COMMON  /RING/  NIBV,RIBB,NIVL,FACTIV,EPIVEL,WFMODL,OPMODL,FGAMMA, 
1  ITERV 

REAL  Q1,W(30),GAMA(30),RS(10),GS(10),WS(10),US(10),DZ{10,36), 

1  R(10,36),DR(10,36),Z(10,36),WIFR(30),U(30) 


COMMON  /HELICOM/Ql,W,U,GAHA,RS,GS,WS,US,DZ,R,DB,Z,WIFR,SHAX 


RS(1)=1.0 

GS(1)=GAMA(I) 

RS(2)=RROOT 

GS(2)=-GS(1) 

DO  IR=1,NIVL+1 

HV=KT( 1 ) *Q1+KT( 2 )*(Q1*IR-Q1 ) 

RT=KT( 4 )+( 1-KT( 4 ) )*( 1 .0/EXP(Ql*IR*KT( 3 ) ) ) 

R(1,IR)=RT 

R(2,IR)=RROOT 

Z(1,IR)=HV 


Z(2,IR)-HV 
END  DO 

DZ  ( 1 ,  NI VL)  =Z  ( 1 ,  NIVL+1 ) -Z  ( 1 ,  NIVL) 
DZ  ( 2 ,  NI VL) -DZ  ( 1  .  NI VL) 

DR( 1 ,NIVL ) -R( 1 ,NIVL+1 ) -R( 1 ,NIVL) 
IF  (WFMODL  .EQ.  1)  THEN 
DR( 1 , NIVL) =1 . 0-R( 1 , NIVL ) 
DR(2,NIVL)-0.0 
DZ(1,NIVL)=0.0 
DZ(2,NIVL)-0.0 
END  IF 
RETLRN 
END 
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APPENDIX  G 
Input  Description 


New  Variables  in  Namelist  NLWAKE 


Variable  Default 
OPMODL  0 


NIVL  4 

WFMODL  2 


FAcrrv  0. 1 

EPIVEL  0.05 

FGAMMA  0.6 

ITERV  200 

NIBV  2 

RIBB(NIBV-2) 


Description 

Inflow  model 

0  vortex  lattice  (existing  CAMRAD  model) 

1  vortex  line 

2  vortex  ring 

Number  of  vortex  levels  in  intermediate  wake;  maximum  36 
Far  wake  model  in  prescribed  wake  vortex  line  and  ring  models 
0  no  far  wake 

1  concentrated 

2  distributed  sheet 

Factor  introducing  lag  in  induced  velocity  iteration 

Tolerance  for  induced  velocity 

Roll-up  weighting  factor 

Maximum  number  of  induced  velocity  iterations 

Number  of  rolled-up  vortices;  minimum  2,  maximum  10 

Inboard  edge  of  rolled-up  vortices  from  root  to  tip,  excluding 

root  and  tip;  must  be  between  root  and  0.9  (default;  evenly 

distributed) 


Comments  on  Important 

Namelist 

Variable($) 

NLTRIM 

MPSl,  MPSIR 

NLTRIM 

DEBUG(14),(24) 

NLTRIM 

LEVEL!  1) 

NLTRIM 

ITERU 

NLTRIM 

ITERR 

NLTRIM 

ITERF 

NLWAKE 

OPRWG 

NLWAKE 

CORE(l),(2),  DBV 

NLWAKE 

FWGT(1),(2) 

Existing  CAMRAD  Variables 

Comment 

Because  solutions  are  independent  of  azimuth,  can  be 
as  low  as  number  of  blades 
Additional  debug  information  for  new  models 
Prescribed  wake  Free  wake 

1  2 

1  Oor  1 

1  Oorl 

N/A  1 

Defines  prescribed  wake  geometry,  as  with  existing 
vortex  lattice  model  in  CAMRAD 
As  with  existing  vortex  lattice  model  in  CAMRAD 
For  Kocurek  &  Tangier  (8  <  OPRWG  <  1 1 )  and 
Landgrcbe  (8  <  OPRWG  S  1 1)  models,  factors  for 
vortex  setding  rates  KT(1)  and  KT(2);  equal  1  for 
unchanged  experimental  rates 
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APPENDIX  H 
Test  Case  Command  File 


$ASSIGN  [USERNAME. CAMRAD.AIRFOIL]NACA0012. TAB  AFTABLEl 
$ASSIGN  [USERNAME. CAMRAD.INPUT1H34. DAT  INPUTFILE 
$DEFINE/USER_MODE  SYS$OUTPUT  [ USERNAME. CAMRAD]H34 .OUT 
SRUN  [ USERNAME. CAMRADJCAMRAD 
&NLCASE  NCASES=1 , BLKDAT=0 , NFRS=-1 , NFEIG=-1 , NFAF2=4 1 , &  END 
&NLTRIM  MPSI=4,MPSIR=4,DEBUG(24)=1,DEBUG(14)=1, 

VKTS=0 . 0 , FACTOR=0 . 5 , ITERC=40 , ITERM=4  0 , FACTM-0 . 5 , MTRIM=80 , 
LEVEL(  1 ) =1 ,  ITERR=1 ,  ITERF=0  ,  DELTA=0 . 5 ,  ITERU=1 , 

OPDENS=3 , TEMP=56 . 0 , DENSE=0 . 00232 , 

OPTRIM=ll , EPTRIM=0 . 001 , CTTRIM-0 . 0817 , MHARM=0 , NROTOR=1 , 
COLL=9.1, 

DOF=54*0,DOFT=8*0, 

ire>RNTP=0 ,  NPRNT1>0 ,  NPRNTT=1 ,  NPRNTI=1 , 

OPREAD=l, 1,2*0, 1,1, 

&END 

&NLRTR  VTIPN=621.6,BT1P=1.00,KHLMDA=1.15, 

KFLMnA=l .  0 ,  FXLMnA=l .  0 ,  FMLMDA-=0 . ,  FACTOU»0 . 5 ,  OPCOMP=l , 

&END 

&NLWAKE  KFW=120,OPHW=0,OPRMG=8,FACIWN=0.1, 

WKMODL(  2  )=0  ,WKMODL(  4  )="0,WKMODL(  6  )=0,WKMODL(  8  )=0  ,WKMODL(  12  )=0 , 
CORE  ( 1 ) =0 . 0 5 , CORE ( 2 ) =0 . 05 , DBV=-1 . , 

FWGT( 1 ) =1 . 0 , FW3T( 2 )=1 . 0, 

EPIVEL=0 . 02 , WFMODL=2 , OPMODL=l , NIBV=2 , FGAMMA-0 . 7 , 

FACTIV-0 . 1 , ITERV«200 , 

&ENO 

&NLBODY  CONFIG=0,ASHAFT(1)«0, 

&END 

&NLLOAD  MALQAD=1,NWKGMP=4*0,HWKGMP=0, 

NPLOT=75*0, 

SEND 
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APPENDIX  J 
Test  Case  Output  File 
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TRI.M  ANALYSIS 

WIND  TUNNEL  TRIM,  OPTRIM 

SINGLE  ROTOR  CONFIGURATION 


OPERATING  CONDITONS 
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NUMBER  OF  HARMONICS  FOR  ROTOR  =  0  lROTOR-1'  0  tR0TOR-2» 

NUMBER  OF  HARMONICS  FOR  AIRFRAME  =  0  (ROTOR-1 ‘  0  (ROTOR-2) 
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